VIBRATION ANALYSIS OF PIPING SYSTEMS 



GEORGE E. FINK 




ubrary 

U.S. NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 



* 3 * 




VIBRATION ANALYSIS OF PIPING SYSTEMS 



* * 



George E» Fink 





1 



I 



VIBRATION ANALISI3 OF 
PIPING SYSTEMS 



by 

George E^ Fink 

Lieutenant, Unltedi^States Navy 



Submitted in partial fulfillment of 
the requirements for the degree of 

i^LASTER OF SCIENCE 

IN 

MECHANICAL ENGINEERING 



United States Naval Postgraduate School 
Monterey, California 

1964 



[fppxsY 

!/‘S. NAVAL P -T. .''aD'J'" fOOCS 

Me ; UE,, C,\Li,., „MA 

VIBRATION ANiiLiSIS GF 
PIPING 3Y3TENS 

toy 

George Eo Fink 



This work is accepted as fulfilling 
the thesis requirements for the degree of 



MASTER OF SCIENCE 
IN 

MECHANICAL ENGINEERING 
from the 

United States Postgraduate School 



aBSTHaCT 



A digital computer program is presented capable 
of determining the in- and/or out-of^plane vibration 
frequencies of a single plane piping systemo The 
analysis is made utilizing the method of transfer matrices, 
a discussion of which is included » The program accepts 
branched systems with or without intermediate supports » 

The effects of both shear deflection and rotational 
inertia are considered although either or both may be 
neglected at the user’s discretiono Distributed mass is 
generally treated as such, although provision is made 
for alternate treatment by lumping masso 

Several typical olping systems are analyzed using 
these various program options, and the results obtained 
are compared* 
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CHAPrER I 



INTRODUCTION 
1 0 1 General Remarks 

The dynamic behavior of piping systems, particularly 
their response to periodic and/or transient excitation, 
is a subject of considerable interest and growing 
pertinence In the design of such systems c The reasons 
for this are readily dlscernable when the unprecedented 
reliability demanded of nuclear power plant piping, or 
the consequence of a piping system failure in an orbiting 
space vehicle, are considered « An increasing number of 
piping systems are being designed, of which, like those 
mentioned above, a long service life is demanded, and 
to which access Is limited or completely denied once the 
system Is placed in service » Clearly, no factor which 
might contribute to the failure of such a system can be 
Ignored* Fatigue brought about by periodic excitation 
at or near a system resonance Is one such factor* Cata-= 
strophlc failure caused by a periodic disturbance of 
short duration and high amplitude, such as an earthquake. 
Is another-* 

With regard to the latter case, it has come to the 
writer’s attention since beginning this work, that 
Japanese engineers are giving considerable and close 
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attention to this particular problenio Their first step 
in approaching it is to determine the natural frequencies 
of the system of Interest o 

Fortunately, or perhaps fostered by necessity, 
development of the “tools'* needed to carry out accurate, 
reliable, and practicable vibration analyses of complex 
systems, namely analysis by matrix methods with the aid 
of high speed digital computers, has kept pace with the 
increasing demand for their employment o 

1,2 Scope of Work Presented 

A method for determining, analytically, the in- and 
out“Of“plane undamped natural frequencies of planar 
piping systems has been developed ana is presented 
herein. It employs the method of "Transfer Matrices"^ 
to construct an accurate mathematical system model and 
determine system natural frequencies, and a digital 
computer to carry out the Inherent voluminous numerical 
computations, VIPIPE, the digital computer program pre- 
sented in appendices A through C constitutes the embodi- 
ment of the method, 

^The theory and method of Transfer Matrices is dis- 
cussed in sufficient depth for the purposes of this 
report in chapter 2, A more comprehensive and detailed 
treatment of the method may be found in ref, 1, 
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This method is a first but important step toward 
the ultimate objective of analysis of three dimensional 
systems of unrestricted topology o 

Attention has been confined to the two dimensional 
or planar case because the theory underlying the treat- 
ment of three dimensional systems using Transfer Matrices 
is not presently available In the literature, and because 
It Is felt that the need for, and utility of, a treat- 
ment of the two dimensional case fully warrants the 
development o 

Although in Its present state of development, the 
theory of Transfer Matrices permits the consideration 
of any planar piping configuration, certain restrictions 
are imposed on the configurations which may be treated 
by the method herein presented, by the limitations of the 
program VIPIPEo These limitations are discussed in detail 
in section lo3‘> 

The philosophy of VIPIPE has sufficient generality 
to make the program readily adaptable to consideration 
of the three dimensional case as soon as the underlying 
theory is developed o 

1 o3 Assumptions and Limitations 

The mathematical models used to describe the 
elastic properties of a system and its components are 
faithful within the framework of the following con- 
siderations s 



3 



a .0 The material of which each component Is made 
is considered to be homogeneous. Isotropic, and 
linearly elastic o 

bo The system is conservative (no damping) o 
Co Distributed mass is treated as such except in 
curved sections for which a lumped mass model is 
usedo For curved sections, sufficient sub» 
sectioning is carried out so as to reduce error 
to an economically^acceptable level o 
do Spring hangers are considered to consist of 
from one to three mutually perpendicular linear 
springs, and from one to three torsional springs 
which act about mutually perpendicular axeso When 
a spring hanger consists of both linear and torsional 
springs, the linear springs are normally considered 
to act along the same mutually perpendicular axes 
that the torsional spring act about o A more gen- 
eral case may be treated with relative ease and 
is a subject of discussion in appendix B, section Bo96« 
Although the theory would permit determining any 
number of natural frequencies for any planar system 
having any number of components, branches, and hangers, 
program VIPIPE, described herein, is limited to the 
treatment of systems having a maximum of fifty components. 



^Appendix B, section B<,93 



twenty branches, and twenty hangers. Further, only 
single filament branches from a main member may be 
treated, i«eo, all branches must emanate from the same 
continuous member, and a branch may not emanate from 
a branch o 
1,4 Notation 



matrix 

x,y,z cartesian right-handed coordinate system 

u,v,w displacements in the x,y,z directions 

respectively 

rotations about the x,y,z axes respectively 
d displacement vector 

E Young’s modulus 

EJ bending stiffness 

G shear modulus 

GJ^ torsional stiffness 

G coordinate transformation matrix 

I unit matrix 

polar second moment of area 
L length 

M bending moment 

N normal force 

0 null matrix 

P internal force vector 

R^^R2 square submatrices 

S spring matrix 
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T torque 

U transfer matrix 

transfer matrix element 

V shear force 

V unknown vector 

Z state vector 

L,R subscripts, left and right side of com- 

ponents respectively 
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CHAPTER II 



TRANSFER 14ATRICES 

2.1 The State Vector and Transfer Matrix 
a. State Vector The state vector at a point 1 of 
an elastic system Is a column matrix having sufficient 
elements to completely describe the Instantaneous dis- 
placement and/or change in orientation of that point 
from its quiescent position, as well as the force and 
moment system existing in the member at the same point 
at the same instant, i.e., those forces and moments 
which, if a cutting plane were passed through the point 
at the Instant of Interest, would have to be applied to 
each cut face to prevent relative motion between them. 

For perfectly general structural systems the 
notion of a state vector is a rather sophisticated 
concept. In any properly constituted state vector 
however the number of elements pertaining to forces and 
moments is equal to the number pertaining to displacement. 

In the case of a planar piping system, the state 
vector has six elements, three pertaining to displace- 
ment, and three pertaining to force. The state vectors 
for the in- and out-of-plane cases ^respectively , as used 

^Henceforth the In-plane case will be referred to 
as Case 1, the out-of-plane case as Case 2. 
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in VIPIPE, were chosen to be 



^lp= 



’U “ 






V 




T 






-w 




^oop- 




Vy 






_N_ 




_Vz. 



(2ol»l) 



The elements of these state vectors are the magnitudes 
of vector quantities of the same name which derive their 
orientations from those of the x,y,z axes of a local 
right-handed cartesian coordinate system (fig« 2ol-l), 
associated with the particular point of Interest » 





Flgc 2ol-l Coordinate system and vectorial representa- 
tion of state vector elements o (Interpret double shaft 
arrows by the right-hand ruleo) 
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These local coordinate systems are constructed 
such that the plane of the quiescent system and the 
xy plane of the coordinate system coincide o The local 
X axis Is coincident with the local tangent to the 
quiescent system centroidal axiSo No reference Is made 
to a common fixed coordinate systemo In establishing the 
coordinate system for the purpose of employing VIPIPE, 
the longest continuous run of pipe, or the one from 
which all branches emanate If they are not one and the 
same, Is designated as the main member, one end of which 
is chosen as the starting point » This end is called 
the left end and its state vector is symbolized o 

The positive x direction is taken as away from the "left” 
end and toward the "right" end along the centroidal 
axis of the main membero The positive z direction is 
perpendicular to the quiescent plane of the system and 
invariant, and this direction is called dowmxardo The 
y axis Is chosen so as to complete a right-handed car- 
tesian coordinate system^ 

b« Transfer Matrix Let us consider the relationship 
between the state vectors at the left and right ends of 
a massless straight section of plpeo We begin by con- 
sidering the situation occurring in Case 2 In which 
both bending and twisting (about the centroidal axis) 
occur. The relationships between the left and right 
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end displacements are given by 



The force relationships are given by§ 
Tr = Tl 

Mr = Ml + v^L 

Vb = Vl 



(2ol-2) 

(2ol-3) 

(2ol=if) 



(2cl«5) 

(2ol=6) 

(2oU7) 



or, in matrix notation 



' 8 ~ 
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0 
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0 


0 
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L 
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The six-by-six matrix relating the state vectors of 
the left and right end of the pipe is termed a transfer 
matriXo It is the transfer matrix of the elastomechical 
system constituted by the pipe and its boundary condi- 
tions o Similar six-by-six transfer matrices may be 
developed for point masses, curved sections of pipe, 
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and elbov;s, etc<, Many of these are contained in the 

1 

transfer matrix catalog Included In ref a 1, and are 
used in VIPIPEo 

Since the application of Interest is to the deter- 
mination of the natural frequencies of piping systems, 
it would be instructive at this point to illustrate the 
origin and nature of frequency dependence In transfer 
matrices,, To that end, let us construct the transfer 
matrix of a point mass, again considering Case 2o To 
begin, note that throughout the application considered 
herein we are concerned with the Isochronous vibration 
of a system at a single frequency C0„ The tacit assiunp- 
tlon has been made that all quantities varying with 
time, l«e», all of the elements of the state vectors, 
do so with the form; constant* coscut, and have exhibited 
only the constant or amplitude, while suppressing the 
variable factor o 

The transfer matrix of a point mass may be constructed 
by noting that deflection, twist, slope, and moment are 
constant and continuous across the zero length section 
which a point mass constitutes « The vibrating mass does 
Introduce a discontinuity in shear however, which may be 
determined by considering the free body diagram of 

^See Appe Fo 
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Fig. 2,1-2, and the following relationships; 



w(t) = w»coso;)t 



Z 2 

<3 w(t) = “WCOCOStjt 



(2.1-9) 



From F=ma: 



Vg-fiCoscJt - V2LC0S0Jt=m§_^^ ^ 

dt^ 



( 2 . 1 - 10 ) 



and substituting (2.1-9) into (2.1-10) we have; 



VzH 



= V^T + mo^(-w) 



( 2 . 1 - 11 ) 



and are ready to construct the transfer matrix Up of 
a point mass. Hence 



( 2 . 1 - 12 ) 






Q 



''zL 



'zTl 



2 

"mo:) w 



Fig. 2.1-2 Free body diagram of a point mass 
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Consider now an elastomechanical system, such as 
that shown in Fig. 2,1-3, made up of many components 
connected end-to-end, each of which, as should now be 
apparent, may be characterized by a unique transfer 
matrix. Subscripting the 




Fig, 2,1-3 An elastomechanical system 

transfer matrices to reflect the subscripts assigned 
to the state vectors which describe the force systems 
and displacements at the ends of each component, we may 
write the following matrix equations; 

^2=Ui, 2*^1> 23=U2^3*225,.o^Zn=Un„i,n'^n-l (2,1-13) 

Note that successive state vectors may be eliminated, 
by substituting for each its matrix product equivalent, 
until finally the matrix relationship 

° remains. The matrix product, 

^n-l,n* • «^3,2*^2,1 » component transfer matrices 

is itself an m-by-m matrix which relates the final to 
the initial state vector of the system. It is, in fact 
the transfer matrix of the system . 
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2,2 Elastic Interaedlate Conditions and Branched Systems 



Piping systems frequently include branches, and 
are supported intermittently along their length by hangers. 
These branches and hangers usually have a very signifi- 
cant effect on the dynamic behavior of a system. Their 
effects may be considered within the framework of the 
transfer matrix method, 

a. Hangers It is clear that the effect of a hanger on 
a system is to introduce discontinuities in the internal 
forces in the system at the point of attachment. 

Consider a spring hanger having the torsional spring 
constants CTX about the x axis, and CTY about the y axis, 
and the linear spring constant CLX in the x direction. 

For a Case 2 situation, the equations relating the state 
vector elements on either side of the point of attachment 
are: 

Tr = Tl + CTX»S 

-Wr =-wl 

Mr = Ml + CTY»Y 

Vr = Vl + CLX-W 
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or In matrix notation; 



's' 




1 


0 


0 


0 


0 


0 




"s 


T 




CTX 


1 


0 


0 


0 


0 




T 


-w 




0 


0 


1 


0 


0 


0 




-w 


¥ 


= 


0 


0 


0 


1 


0 


0 




¥ 






0 


0 


0 


CTY 


1 


0 




% 




fi 


0 


0 


^CLX 


0 


0 


1 




V2 



It should be apparent at this point that a rigid 
intermediate support, ioSo, one having one or more 
infinitely stiff springs, may not be considered using 
this method o Rigid supports may be approximated however 
by using very large values for appropriate spring con- 
stantso A value of 100,000 IbSo/ino or ln»=lbo/ rad» 
as applicable would be an adequate approximation in most 
instances, and it does not introduce computational diffi- 
culties , 

bo Branches The Joining of a branch to a system 
has an effect similar to that of a spring hanger on the 
force system in the main member at the point of Joining «. 

In fact, it "looks" exactly like a spring hanger to the 
main member except, as vjill shortly be seen, its parameters 
depend on frequency* Discontinuities are Introduced in 
the forces whereas the main member deflections on each 
side of the Joint are the same, and Identical to those 
of the branch at the same point* 
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Consider a Case 1 situation and systeji shown in 



Figo 2o2-l. 




Figo 2o2“l Main system with a branch 
In the free body diagram, Figo 2o2-2a, of the joint of 
the system of Figo 2o2-l, the force systems acting on 
the Joint are expressed in terms of their own coordinate 
systems 0 ^ By a suitable coordinate transformation it 
is possible to express the force system of the end of 
the branch in terms of the main 




member coordinate system o 
Fig. 2.2-2b. 




Figo 2o2-2b 

Such a situation is shoi-m in 



iThe symbol " refers to the fact that the vector 
elements are expressed in terms of the coordinates of 
the branch 0 
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Equilibrium considerations and the orthogonality of 

the force and shear vectors permit us to write the 

following equations: 

M = M + M 
RLE 

V = V + V 
RLE 

N = N + N 
RLE 



or in matrix form 











M 


V 


+ 


V 


= 


V 






N 




N 



L B R 

Bearing in mind the continuity of displacements 

at a joint, it is clear that the main member state 

vector Zg can be obtained by simply adding the force 

components of the branch state vector Zg, expressed 

in the coordinate system of the main member, to the 

force components of the main member state vector Z^<, 

What follows is the derivation of a suitable expression 

for the force components of the branch state vector 

The main member state vector is obtained easily 

from the state vector Z„ and the transfer matrices of 
A 

the Intervening main member components o 

At point C (Flgo 2 o 2-1) there are three zero and 
three non zero elements in the state vector 

^See Appendix B, section Bo 92 for treatment of end 
conditions which do not meet the requirement for three 
zero and three non zero elements o 
17 
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The coluom matrix of the non zero components is 
symbolized Vo after relating the branch state vector 
at B to that at C by the transfer matrix of the branch 
(obtained in the usual manner from the transfer 
matrices of its components), and eliminating those 
columns of the branch transfer matrix which would 
otherx-flse be multiplied by the zeros of Zq, we obtain 
the matrix equation 







fR' 






1 


"1" 


— 


— 


Pb 







wherein the state vector Zg has been reordered and par- 
titioned such that ^denotes a sub vector containing 
the elements of Zg pertaining to deflection, and ^a 
sub vector containing those pertaining to force o Sxpan 
Sion of equation 2o2-l gives 

dQ = Ri o V (2o2-2) 

and 

Pb = R 2 » (2c2-3) 

By substituting “ <lg = v In equation 2o2-3 we 
obtain 

h = h’ ” 5b (2.2-4) 

relating the displacements at B to the forces in the 
branch at Bo The matrix product is called a 

spring matrix because of the function that it serves 0 
18 
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iilthough any consistant system ..ill do, it is con- 
venient here to establish a particular method for measuring 
the angle at which a branch Joins the main membero The 
one defined is consistant with that used in VIPIPEo 
Accordingly, the angle at which a branch Joins the main 
member is measured counter-clockwise from the x axis of 
the main member as it is oriented immediately before 
the intersection, to the x axis of the branch as it is 
oriented at the point of intersectlono The angle 0 in 
Figo 2o2-l is drawn in conformance with this deflnltlono 
From the matrix relationships 

dg = Gi“ dg 

and 

I’B = <52* 

where for a planar piping system 



COS0 


-sin0 






1 


0 


0 


sin0 


COS0 


0 


and Gg = 


0 


COS0 


-sln0 


0 


0 


rj 




0 


sln0 


COS0 



the spring matrix (symbolized 3) < 

S = Gg- 

is obtained which relates the force and displacement 
elements of the branch transfer matrix Zg (expressed in 
terms of the coordinates of the main member) o Thus 

Pb = 3" dB 
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It is now a simple matter to construct the transfer 
matrix relating the main member state vectors on either 



side of the joint o Recall from equation 2»2-0 that 



which has just been derived o 

Thus a branch is treated as a point transfer matrix 
in constructing the transfer matrix of the system® 

2.3 Mode Frequencies and the Frequency Determinant 

A principal application of the method of transfer 
matrices, and the one x\rith which we are concerned, is 
to the vibration analysis of elastomechanical systems 
for the purpose of obtaining their natural frequencies® 
When transfer matrices are constructed for systems 
having mass, certain of their elements are frequency 
dependent (recall the transfer matrix of a point mass 
derived earlier)® Proper application of the system 
boundary conditions to these transfer matrices yields 
a matrix of order n/2 (from a transfer matrix of order n) 
which can be used to obtain the natural frequencies of a 






since 



^R ■ , 




I I 0 fd" or = Up. 



where we use the symbol to denote the point matrix 
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system. Its determinant is termed the frequency determi- 
nant of the system. The concept is best described by 
an example. 



Consider a case 2 situation wherein a system has 
a system transfer matrix U 3 and boundary conditions 
embodied in its left and right end state vectors and 
Zjj. Assume that the left end is fixed and that the right 
end is free. The matrix equation relating the left and 
right end state vectors Is 



's' 




‘ui,i 0 


. 




""1,6 




O' 


0 






. 




0 




T 


-w 




. 


, 


0 






0 


0 




0 


, 








0 


.oj 


R 


_U5^1 . 


, 




^^6,6 




h 



from vrhich may be constructed the matrix relationship 



'o 




^2,2 


""2,5 


^2,6 




T 


0 


= 


^5,2 


^5,5 


^5,6 


• 




0 


R 


^6,2 


^6,5 


^6,6 







Clearly, a non-trivlal solution requires that the third 
order submatrix of equation 2,3“1 must vanish. The 
determinant of the submatrix is the system frequency 
determinant mentioned above. The frequencies for which 
its value is zero are the natural frequencies of the 
system described by the transfer matrix U 3 and the 
boundary condition state vectors, Z^ and Zjj, from which 
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it was derived 



2 o4 The Dlpqltal Computer and the Method of Transfer Matrices 
Even for a relatively simple system, the elements 
of the frequency determinant become unwleldly polynomials 
in radian frequency o Determination of the roots of the 
equation resulting from an expansion of the frequency 
determinant is usually quite impractical and recourse 
must be made to a numerical solution » In a numerical 
solution the frequency determinant is constructed and 
evaluated for selected frequencies over a range of 



interest 0 Its value is plotted as a function of fre- 
quency and the natural frequencies of the system which 
it represents are the zeros of the plot (Figo 2<,3-l)o 




Use of a digital computer permits rapid and accu- 
rate analysis of complex systems « Program VIPIPE 
presented in Appendices A through C, for instance, will 
accept for analysis planar piping systems having up to 
fifty components Including as many as twenty branches 
and spring hangers. It treats distributed mass as such 
and considers the effects of shear deflection and rota- 
tional inertia. In a typical analysis, the first five 
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natural frequencies were found for a system having 
thirteen components, including two spring hangers and 
two branches, at the cost of five minutes of machine 
time on a GDC l604 computer c 

It is throup-h the use of high speed digital 
computers, that the full potential of the method of 
transfer matrices in this application as well as in 
others^may be realized. 



^See ref. 1 for a comprehensive treatment of 
applications of the method of transfer matrices to 
other problems in elastomechanlcs . 
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CHAPTER III 



DISCUSSION 



3.1 Accuracy of Method 

Confidence in the accuracy of the method was 
established by using it to determine the natural 
frequencies of several systems for which these fre- 
quencies could be obtained by other means, loeo, by 
means of either a closed form analytical solution or 
by recourse to the literature The average dif- 

ference between system natural frequencies obtained 
using VIPIPE and those obtained by analytical means 
was found to be less than .0064^, while the average 
difference between system natural frequencies obtained 
using VIPIPE and those obtained from the literature 
was found to be less than 

Details of the systems analyzed and the results 
obtained may be found in Appendix E, section Eo21. The 
analyzed systems consisted of single spans of straight 
pipe having various combinations of boundary conditions o 
The published papers of E» Voltera [ 4 ] give 
methods for determining the natural frequencies of 
curved filaments, but the methods are sufficiently 
Involved as to make their application in obtaining 
some comparison frequencies a major undertaking., 

^Numbers in brackets pertain to bibliography. 

24 



Time precluded assumption of the task. The literature 
[2,J does provide means for computing the fundamental 
frequency of partial rings with fixed ends » The dif- 
ference between the fundamental frequencies determined 
using VIPIPfi and those computed using the method of 
the literature, was found to be less than, 3;^ for in- 
plane, and 2,3^ for out-of-plane vibrations Details 
of the system analyzed and the results obtained may be 
found in Appendix E, section Eo24o 
3.2 Integrity of Solution 

Confidence in the integrity of the method was 
established by conduct lng”mirror image" analyses of 
almost every configuration considered. Construction 
of the transfer matrix of a system is directional in 
nature, in that an end of the system is chosen as a 
starting point, and the transfer matrix of the system 
is developed by premultiplying the Incomplete system 
transfer matrix by the transfer matrix of each system 
component in the order that each component is encoun- 
tered in tracing the system from that starting point. 

It follows therefore that an error in the structure of 
any of the component transfer matrices, or in the method 
or order of their combination in constructing the 
system transfer matrix, would reveal itself in the form 

^See footnote to section E,24- of Appendix E, 
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of inconsistent results if the system were analyzed 
twice, beginning the first time from one end, the second 
time from the other endo Such a dual analysis is 
referred to herein as a “mirror image” checko 

In general, the disagreement between natural 
frequencies obtained for a system and for its mirror 
image was found to be less than o 0003 ^<> 

An added check on solution integrity was obtained 
by analyzing a single branch system treating first 
one member then another as the branch o The average 
difference here was on the order of .02?^ for in-plane 
vibration, and .00?^ for out^of^plane vibrationo 

Details of the systems analyzed and of the results 
obtained may be found in Appendix E, sections Eo31 and E,32. 
3.3 Computer Imposed Limitations 

The number of natural frequencies that may be found 
for a given system using VIPIPE, is a function of the 
floating point significant figure capacity of the com- 
puter usedo The significant figure capacity required 
to solve for a given number of natural frequencies of 
any system is a function of the characteristics the 
system; and for all systems, the significant figure 
capacity required increases with the number of the mode 
sought. Thus, for example, for a given system, the 
number of significant figures required in order to 
determine the third mode frequency is greater than 
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that required to determine the first mode frequency, 
and vjhile the first ten natural frequencies of one 
system may be obtained without difficulty, the signifi- 
cant figure capacity of the computer may be inadequate 
to determine the fourth mode frequency of another* 
Program VIPIPE is designed to terminate a solution 
when the number of significant figures required for its 
continuance exceeds that available in the machine* 

The use of multiple precision arithmetic would 
extend VIPIPE'S utility considerably* 
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APPENDIX A 



DESCRIPTION OF PROGRaM 
Aol General Remarks 

Program VIPIPE is a CDC FORTRAN 60 language 
digital computer program designed for use in the vibra- 
tion analysis of planar piping systems c In and out-of- 
plane undamped natural frequencies may be found through 
its use. The analysis is performed by an application 
of the method of transfer matrices wherein each component 
of a system is characterized by a matrix called its 
transfer matrix. The entire system is characterized 
by the system transfer matrix obtained through an 
ordered multiplication of the component transfer 
matrices. Considerable flexibility with regard to 
the mathematical model employed to describe the system, 
and the nature of simplifying assumptions accepted in 
its vibration analysis, is provided in the program. 

For Instance, although the method of solution yielding 
the greatest accuracy coupled with least machine time 
combines a distributed mass treatment of straight 
sections with a lumped mass treatment of curved sections, 
and considers the effects of both shear deflection and 
rotational Inertia, the user may elect to use a lumped 
mass treatment of all sections and to consider or neglect, 
singly or together, the effects of shear deflection and 
rotational inertia. 
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Problem solution proceeds by an Iterative process 
with radian frequency the independent variable « The 
sign of the frequency determinant^ is used to control 
the search for and convergence to the system natural 
frequencies. The frequency Increment is computed 
internally, as is the solution acceptability criterion. 
Both are accessible, however, and may be controlled by 

p 

the user. Subsectioning of components treated by lumped 
mass methods is carried out Internally as a function of 
component L/D (length-to-diameter ) ratio and the highest 
mode sought. 

The program has the capability to analyze piping 
systems comprised of a maximum of fifty components, and 
having as many as twenty branches and twenty hangers , 

A component, for program purposes, is defined 
as a curved section, i,e,, a section of constant curva- 
ture, a straight section, a spring hanger, a flange or 
a coupling,^ Although straight and curved sections, 
flanges and couplings must each be of constant outside 
diameter and wall thickness, these dimensions may vary 
from component to component. Similarly, intensive 
properties, i,e,, shear and elastic moduli, and density, 
may not vary within a component, but may vary from com- 
ponent to component, 

1 

See Chap, 2, 2,3 of the main body, 

2see App, B, section B,94, 

^See App, B, section B,95« 
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A shear area form factor of one half is built into 
the programo It may be altered at the user’s discretion 
by altering one cardc^ 

Five common end or boundary conditions, for the 
ends of the main member and the remote ends of any 
branches, may be specified by a simple coded input* 

Unusual boundary conditions may be considered by recourse 
2 

to a more complicated input® 

Systems having rigid Intermediate supports and/or 
non linear elements may not be analyzed by the use of 
this program® The number of natural frequencies found 
for a given system is a function of the significant 
figure capacity of the computer employed since the num- 
ber of significant figures required in a solution is a 
function of the characteristics of the system under 
analysis, and for any system Increases with increasing 
mode number® 

A® 2 Program Structure 

The program consists of a six section main body and 
twenty-four subroutines® If sine, cosine, slnh, and 
cosh function routines are not available in the FORTRAN 
library they must be provided by other means since, 
although required by the program, they are not explicitly 
contained therein® The function of each section of the 
main body and of each subroutine is given below. 

^See App® B, section B®91« 

^See App® B, section B®64® 
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Ao21 Main Body Sections and their Functions 



IKPUT Control read-ln of data, allocate storage 

locations for arrays and subscripted vari- 
ables 0 



IKVaRIaNTS Compute, and assign appropriate subscripts 
to, such system component characteristics 
as moment of Inertia, and radius of gyra- 
tion. 



CONTROL 1 Construct transfer matrix of the system 
for the in-planelcase by calling for 
various subroutines in a sequence governed 
by the system geometry, the mathematical 
model specified, and the desired assumptions, 

CONTROL 2 CONTROL 2 performs the same function as 
CONTROL 1 but for the out-of-plane^case. 



ITERATION Apply boundary conditions to the system 

transfer matrix constructed by CONTROL lor 2, 
so as to obtain the system frequency 
determinant. Evaluate the determinant 
and utilize its slgn^ln controlling the 
Iterative search for mode frequencies. 
Initiate output when the specified fre- 
quency range has been fully explored, the 
required number of modes found, or the 
significant figure limit of the comouter 
reached. 



OUTPUT Control output format. 



A, 22 Subroutines and their Functions 



SUBSEC Subdivide components to be treated by 

lumping mass on the basis of L/D ratio 
(or ratio of radius of curvature to diam- 
eter ^and Included angle of arc, as 
applicable , ) 

DISTM Construct transfer matrix of straight' 

sections having distributed mass for Case 1, 



^Henceforth Case lo 
^Henceforth Case 2, 
pSee App, B, section B,9A, 

^Henceforth RHO/D where RHO is the radius of curvature 
of the component. 



DI3TM0 


Construct transfer matrix of straight 
sections having distributed mass for 
Case 2o 


3 FI ELD 


Construct the transfer matrix of straight 
massless subsections for Case 


3FIEL0 


Construct the transfer matrix of straight 
massless subsections for Case 2o 


CFIELD 


Construct the transfer matrix of curved 
massless subsections for Case lo 


CFIELO 


Construct the transfer matrix of massless 
curved subsections for Case 2, 


POINT 


Construct the transfer matrix of a point 
mass for Case lo 


POINO 


Construct the transfer matrix of a point 
mass for Case 2 <> 


MATMUL 


Construct the transfer matrix of those 
components treated by lumping mass through 
the ordered matrix multiplication of the 
transfer matrices of massless subsections 
and point masses o 


RIGID 


Construct the transfer matrix of a rigid 
straight section having distributed mass 
for Case 1« 


RIGIO 


Construct the transfer matrix of a rigid 
straight section having distributed mass 
for Case 2» 


3TIFC0 


Construct the transfer matrix of a massless 
rigid corner (essentially a transformation 
of coordinates at an elbow) for Case 1« 


STIFOO 


Construct the transfer matrix of a massless 
rigid corner (essentially a transformation 
of coordinates at an elbow) for Case 2» 


HANGER 


Construct the transfer matrix of a spring 
hanger for Case 1<, 


HANGEO 


Construct the transfer matrix of a spring 
hanger for Case 2« 


BRANCH 


Construct a transfer matrix from the trans- 
fer matrix of a branch and its remote end 





boundary conditions, i^hich permits a 
branch to be treated as a component of 
the main member ^ for Case lo 


BRaNCO 


Construct a transfer matrix from the 
transfer matrix of a branch and its 
remote end boundary conditions, which 
permits a branch to be treated as a com- 
ponent of the main member for Case 2 o 


STAVEC 


Construct boundary condition vectors 
from coded input for Case lo 


STiiVEO 


Construct boundary condition vectors 
from coded input for Case 2 » 


INVERT 


Invert a matrix. 


FINHAT 


Construct the transfer matrix of the 
system by successive premultiplication 
of the partial system transfer matrix by 
the transfer matrix of each system com- 
ponent as soon as each of the latter is 
constructed 0 


FINBRA 


Construct the transfer matrix of a branch 
by successive premultiplication of the 
partial branch transfer matrix by the 
transfer matrix of each branch component 
as soon as each of the latter is constructed 


DELMA 


Construct the frequency determinant of 
the system by applying the boundary con- 
ditions of the main member to the com- 
pleted system transfer matriXo 



Ao 3 Propiram Nomenclature 
A. 31 Main Body 



A 


Array /matrix 0 Transfer matrix of a 
massless subsection. 


AA 


Discriminant 0 A control in the conver- 
gence sequence. 


AAMU 


In INPUT; density ( lbs./cu , f t , ) of the 
first system component. 



In IKVARIAWTS; average mass per unit 
length of the main member 1- 

aMU In INPUT; component density (lbs,/cUo fto)o 

In remainder of program; component mass 
per unit length o 

aD Average outside diameter of the main 

member 

iiDOi'l Current computed Iteration frequency 

increment o 

AFR Average moment of inertia of the main 

member 0 3 

AE Elastic modulus of the first system 

component o 

aG Shear modulus of the first system com- 

ponent o 

aid Average inside diameter of the main member^ 

AIX Component radius of gyration about the 

X axiSo 

aIY Component radius of gyration about the 

y or z axiSe 

AJ Component moment of inertia about the 

y or z axiSo 

AJT Component polar moment of inertlao 

AL Total length of main member « 

AiiYK Discriminant o Controls choice of mathe- 

matical model, i^eo, lumped or distributed 
mass 0 

B MatriXo Transfer matrix of a point masso 

BB Discriminant o Controls iteration fre- 

quency during convergence « 



^Employed in computing initial frequency increment and 
solution acceptability crlteriono 
2-^Employed in computing the initial frequency Increment 
and solution acceptability crlteriono 
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BDOM 


Initial frequency increment * 


BAriU 


Component density (lbs./cUo ft*) 


3H 


Muinber of last mode found* 


BR 


Initial solution acceptability criterion. 


CC 


Discriminant. Designates final iteration 
during convergence. 


CLX 


Linear spring constant of a hanger in 
the X direction.^ 


CLY 


Linear spring constant of a hanger In 
the y direction.^ 


CLZ 


Linear spring constant of a hanger in 

the z direction. 3 


CR 


Current value of the solution acceptability 
criterion. 


CTX 


Torsional spring constant of a hanger about 
the X axls.^ 


CTY 


Torsional spring constant of a hanger about 
the y axis. 5 


CTZ 


Torsional spring constant of a hanger about 
the z axis.® 


D 


Component outside diameter. 


DD 


Discriminant. Controls read-ln and stor- 
age of the elastic and shear moduli and 
density of components. 


DI 


In INPUT; component wall thickness. 

In remainder of program; component inside 
diameter. 


DT 


Component wall thickness . 


DOM 


Current frequency increment during con- 
vergence. 



^"^No coupling is permitted in action of hanger* For 
treatment of hangers in which coupling occurs see 
App. B, B.9. 



Component elastic modulus o 

Average elastic modulus of the main 
member. 

Shear area form factor. 

First mode frequency of a fictitious 
straight beam having fixed ends, physical 
properties and transverse dimensions 
equal to an average of the characteristics 
of the main member to which it Is equal 
in length. 

Component shear modulus , 

Discriminant, Control In test for approach 
to significant figure limit of computer. 

Reference value of the frequency deter- 
minant (significant figure limit test). 

Number of highest mode sought. 

Subscript , 

Subscript which associates a branch with 
the angle at which it joins the main 
member. 

Discriminant, Controls analysis to pro- 
vide Case 1 and/or Case 2 mode frequencies 
as specified. 

Discriminant. Controls solution sequence 
during Case 1 or 2 analysis. 

Discriminant, Controls output to provide 
a print-out that an analysis has been 
terminated due to significant figure limi- 
tations when appropriate. 

Iteration number. 

Number of current iteration. 

Number of mode currently being sought. 

Number of highest mode found. 

Subscript which associates a branch with 
its remote end boundary condition. 



WBR 


Number of branches in a systera. 


NM 


Discriminant c Controls output* 


NMBR 


Total number of system boundary conditions 


NN3R 


Discriminant* Assigns a component to 
the main member or to a branch, and if 
the latter, indicates its position in 
the branch* 


NS 


Total number of system components* 


Ol-iEGA 


Current iteration frequency* 


OiiEGAC 


Convergence reference frequency* 


OMEGAG 


Upper limit of frequency range of investi- 
gation* 


OWEGAI 


Investigation starting frequency* 


OMEGAH 


Mode frequency* 


P 


Sum of the numerical values of the spring 
constants of a component * 


PDSL 


Value of the frequency determinant used 
as a reference during iteration* 


PHI 


Angle at which a branch Joins the main 
member. 


R 


Reciprocal of component flexural rigidity. 


RHO 


Component radius of curvature* 


RI 


Discriminant* Controls consideration 
of rotational inertia. 


SD 


Discriminant* Controls consideration of 
shear deflection* 


3V 


Coded boundary condition for Case 1. 


SVIN 


Non-coded boundary condition for Case 1* 


SVO 


Coded boundary condition for Case 2* 


3V0P 


Non-coded boundary condition for Case 2. 


TDEL 


Discriminant* Controls significant figure 
cut off* 
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THiiTA 

THETR 

TL 

U 

uu 

V 

w 

z 



Included ansle of arc of curved com- 
ponents in degrees in input, radian 
measure in remainder of program. 

Included angle of arc of curved com- 
ponents in degrees. 

Total length of a component. 



Array name. 

Array name. 

Array name, 
matrix. 

Array name. 



Component transfer matrix. 
System transfer matrix. 
System or branch transfer 

Branch transfer matrix. 



Number of subsections of a component. 



A. 32 Subroutines 



a. SUBSBC 



RHOD 


Absolute value of radius of curvature of 
component . 


RA ■ 


Component L/D ratio. 


RD 


Component RHO/D ratio. 


b. DITii and 


DISTNO 



ara 



Component cross sectional area. 



AREA 

ALl^ 

AL2 

CBS 

CHLl 

CL2 

COO 

COl 

C02 

C03 

BE 

CLl 

S 

SBE 

SHLl 

SL2 

SPT 

W 






Component equivalent shear area. 



Array elements and/or computed constants 
used in evaluating array elements. 
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Component flexural rigidity. 

Component torsional stiffness (DISTFiO). 
and SFIELO 



Component cross sectional area (3FIELD only). 
Component equivalent shear area. 

Length of subsection. 
and CFIELO 



Component cross sectional area. 

Component equivalent shear area (CFIELD only). 



Array elements and/or constants used 
in computing array elements. 



W Component shear stiffness (CFIELO only). 

PHI Included angle of arc of subsection. 

RD Absolute value of radius of curvature 

(CFIELD only). 

W Component longitudinal stiffness (CFIELD 

only ) . 

RHOD Signed component radius of curvature 

(CFIELO only), 

SD Discriminant, Controls consideration 

of shear deflection (CFIELO only), 

V Factor. Imparts sign to certain array 

elements , 

e, POINT. POIKO. PIGID. and PIGIO 

AM Total mass of subsection or component. 

f . STIFCO and STIFOO 



V 



Factor, Imparts sign to certain array 
elements . 
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g. BRANCH and BRiiNCQ 



W 

uv 

SVB 

SR 

Rl, 

R 

S 

Gl, 

U 



Branch transfer matrix (BRANCH) 

Branch transfer matrix (BRilNCO). 

Array name. Remote end boundary condi- 
tion vector (BRANCH) . 

Array name. Remote end boundary condi- 
tion vector (BRANCO), 

R2 Array names. Submatrices of branch 

transfer matrix. 

Array name. Matrix product of R2 " Rl”^ 
in first usage, spring matrix of branch in 
second usage in BRANCH, 

Array name. Spring matrix of branch 
BRANCO . 

G2 Coordinate transformation matrix elements. 

Array name. Transfer matrix which repre- 
sents a branch as a component of the main 
member , 
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APPENDIX B 



INSTRUCTIONS FOR PROGRAM USE 
B. 1 General Remarks 

Appendix B is intended to provide complete instruc- 
tions for using program VIPIPEo To that end, since for 
problem programming purposes, most system components must 
be identified with respect to a coordinate reference, a 
description of the reference system employed is given first 
It is followed by a summary of data required, a description 
of the component identification system, an explanation of 
the nature of a state vector in terms of whose elements 
the system's boundary conditions, and its displacements 
and internal force systems are described, and a description 
of the data input format. The appendix concludes with 
some comments on the options available in the program. 

For the purpose of setting up a problem for solution, 
the longest continuous run of pipe, or the one from which 
all branches emanate if these are not one and the same, 
must be designated as the main member. A system component 
is defined as a section of pipe, straight or of constant 
curvature, an elbow, a spring hanger, a coupllng^or a 
flange. Diameter, wall thickness, elastic and shear moduli 
and density may vary from component to component, but may 
not vary within a component. Since both torsional and 
linear spring rates (for hangers) appear Internally as 

^See App, B, section B.95« 

B=^l 



ele^ients of a matrix used in obtaining matrix products, 
they must have finite values c Thus perfectly rigid 
hangers may not be treated exactly, but they can be 
approximated by using very large spring constants* 

B.2 Coordinate System 

In accordance vjith the procedures given below, a 
local coordinate system is constructed at each point of 
interest in the system. These points are usually the 
ends of each component, the beginning and end of the main 
member, and both ends of each branch. The coordinate 
systems are established with the piping system quiescent. 
Once established, they remain fixed in space. 

In establishing a coordinate system, one end of the 
main member, henceforth called the left end, must be chosen 
as a starting point. 

A right-handed coordinate system is used whose x 
axis, always directed away from the left end, is the local 
tangent to the central axis of the main member. The origin 
of the coordinate system lies on the main member central 
axis, and the xy plane coincides with the quiescent plane 
of the system. The z axis is always directed in the same 
direction, and this direction shall be called downward. 

The direction of the local y axis is chosen so as to com- 
plete the right-handed system. 




Z 

Fig. B.2 Coordinate System 
B“2 






* 





The starting point, for reference purposes, of a 
branch is taken as its far end, ice,, the end farthest 
removed from the intersection of the branch vjith the main 
member o The plane of the branch coincides with the xy plane o 
The positive z axis of the branch and main member are 
parallel and have the same sense. The y axis is oriented 
so as to complete a right-handed system, 

B.3 Branch Joining Angle 

The angle at which a branch Joins the main member is 
measured counter-clockwise in degrees from the x axis of 
the main member as it is oriented immediately before the 
intersection, to the x axis of the branch as it is 
oriented at the point of intersection, 

B.^ Information Required to Program Problem 
B.41 General 



a, number of highest mode sought, 

bo frequency range over which investigation is 

to be performed (initial and final frequencies 
in radians per second), 

c, vjhether in-plane^ out-of-plane or both in- and 
out-of-plane mode frequencies are sought, 

d, whether employment of a distributed or of a 
lumped mass model is desired, 

e, whether it is desired to consider or ignore 
singly or together the effects of shear 
deflection and rotational inertia, 

fo boundary or support conditions at the ends 
of the main member and at the remote end of 
each branch, 

g, number of components in the system, 

h, number of branches in the system. 
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B.42 



All Coaponents ^xcent Hanp-ers 



a, outside dlajieter ( Inches )o 
bo wall thickness ( inches), 

Co density (lbs«/cUo f t>) , 
do shear modulus (psi), 
e» elastic modulus (psi)o 
B.43 StralA’ht Sections 

a, length (inches), 

B.^4 Curbed Sections 

a, signed radius of curvature (negative if 
center of curvature lies on the negative 
y axis, positive if on the positive y axis) 

bo arc central angle (degrees )o 

B 0 ^5 Hangers ^ 



a, linear spring constant In x direction 
(lbs o/ino ) o 

b, linear spring constant in y direction 

(IbSo/ino ) o 

Co linear spring constant In z direction 

(IbSo/ln, ) o 

do torsional spring constant about x axis 
(in«-lbSo/rado ) o 

e» torsional spring constant about y axis 
(in«-lbSo/rado ) , 

fo torsional spring constant about z axis 
( In, “lbs o/rado ) o 

B.5 Component Identification System 

Components are identified by two numbers. The first 
a component sequence number, Is used to order data cards 
containing component extensive and Intensive property data 



ISee iippo B,, section 3,96, Chap, 2, section 2o2a, 



The second, a branch component identification number, is 
read in as data and serves to indicate :;hether or not a 
component is a member of a branch and if so, the position 
it occupies in the branch o 

In addition, all branches, if any, are assigned a 
sequence number which is used to order data and/or data 
cards pertaining to branch boundary conditions and Joining 
angles. 

The methods of assigning these numbers are given 

below, 

B.51 Component Sequence dumber 

All system components bear a component sequence 
number. The number one is assigned to the leftmost com- 
ponent of the main member. Each succeeding component 
encountered in tracing the main member is assigned a 
number one greater than that assigned to its lefthand 
neighbor. If and when a branch is encountered its com- 
ponents are also assigned sequence numbers beginning with 
the component farthest removed from intersection of the 
branch with the main member. That component is assigned 
a number one greater than that assigned to the last com- 
ponent of the main member encountered prior to, and abut- 
ting on, the Junction, When more than one branch Joins 
the main member at the same point, the user is free to 
choose the order in which desired to take up the branches 
for the purpose of assigning sequence numbers to their 



components. Having established the order, the first 
branch is treated as if it xvere the only branch. The 
remaining branches are treated in exactly the same manner 
except that the component of the preceeding branch having 
the highest sequence number, i.e., the one abutting the 
joint, assumes the role performed by the last encountered 
main member component in treating the first branch (see 

Fig. B.51). 




Flgo B51 Component Sequence Numbers 
B«52 Branch Component Identification Number 



Assignment of a component to a branch is indicated 
by the value of a fixed point constant (NNBR) coded as follows 

a. NNBR=1 if component is first member of a 

branch, i.eo, the component 
farthest removed from the inter- 
section. 

b. NNBR=2 if component is neither the first 

nor the last member of the branch. 

c. NNBR=3 if the component is the last 

member of the branch. 
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do NK3R==1 if the branch has only one com- 
ponent o 

e» NNBR=0 all main member components «, 

Bo 53 Branch Sequence Number 

Branches are mimbered sequentially as they are 
encountered In tracing the main member® If tx'jo or more 
branches join the main member at the same point, the branch 
whose components have been assigned the lowest sequence 
numbers is assigned the loiter or lowest branch sequence 
number (see Flgo Bc53)« 




Flgo B.53 Branch component identification 
and branch sequence numbers® 

B 0 6 Boundary Condi t ion Codification 

A boundary condition of the main member, or of a 

branch if any, consists of six parameters which, in the 

program, constitute the elements of a column matrix called 

a state vector® Five common boundary conditions, those of 

a fixed, free, pinned, propped, or a roller supported end 

may be indicated by a code given below which the program 

can Interpret, and from v;hlch it constructs the required 

state vector.. 



o.ol Boundary Condition Code 



Codable boundary conditions are Indicated by the 
below listed values of the floating point constant SVs 
ao SV=le fixed end 
bo SV=2o free end 
Co 3V=3o pinned end^ 
do SV=4o propped end^ 

3 

eo SV=5o roller supported end*^ 

'i/^hen a boundary condition Is other than one of the 
five expressable by this code, the user must construct the 
necessary state vector and cause it to be read in as data» 
This situation will almost always obtain when both Case 1 
and Case 2 mode frequencies are sought for the same system. 
It Is apparent, for Instance, that a constraint which fits 
the definition of a prop when considering Case 1 vibrations, 
may not be descrlbable by any of the codes available when 
Case 2 vibrations are considered. 

In constructing a state vector it is only necessary 
to know whether the value of a constituent parameter is 
zero or non-zero. The code for a parameter having a value 
of zero is Oo, while that for one having a non-zero value 
is 1.0 

The constitution of the state vectors for Cases 1 
and 2 are given below. 
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See App, B, Bo64 Catalog of Coded State Vectors 



Bo 62 In-Plane State Vector (Case 1) 



The vector elements In order are; 
u displacement In the x directiono 
V displacement In the y directiono 

l9' slope, In the xy plane, with respect to the x axis 
moment about the z axis« 

Vy shear In the y directiono 
N normal force In the x directiono 
B«63 Qut-of-Plane State Vector (Case 2) 

^ twist about the x axiSo 
T torque about the x axiSo 

-w (negative of) deflection In the z directiono 
Y slope, in the xz plane, with respect to the x axis 
My moment about the y axis o 
V 2 shear in the z directiono 




Figo B 063 Vectorial representation 

of state vector parameters 



B 0 64 Catalo/a; of Coded State Vectors 



a o___Ca^e I s 







Boundary Conditions 


t 


State 

Vector 

Parameter 


















"propped 


p 


fixed 


free 


pinned 


roller support 


u 


Ou 


lo 


Oo 


lo 


Oo 


V 


0, 


lo 


Oo 


Oo 


lo 


tS- 


Oc 


lo 


lo 


lo 


lo 




Ic 


Oo 


Oo 


Oo 


Oo 




lo 


Oo 


lo 


lo 


Oc 


N 


1. 


Oo 


lo 


Oo 


lo 



bo Case 2o 







Boundary Conditions 


t -■ 


State 


1 










Vector 

Parameter 


1 

fixed 


free 


"'ginned 


'^^^opped 


roller support 


s 


Oo 


lo 


Oo 


lo 


lo 


T 


lo 


Oo 


lo 


Oo 


Oo 


-w 


Oo 


lo 


Oo 


Oo 


lo 


T 


0. 


1 0 


lo 


1, 


1, 


“y 


lo 


Oo 


Oo 


Oo 


Oo 




lo 


Oo 


1 0 


lo 


Oo 
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Problem Statement 



The first data card contains a statement of the 
problem» The card has ten fields under the control of the 
following field specifications: 

312 , 4F3oO , 2P8oO , F2.0 
Field names and the information each field conveys 
are given below followed by an explanation of codes used. 
Field Name and Content 



Ic NS 


number of components o 


2 . NBR 


number of branches o 


3 c lOP 


code calling for in and/or out-of-plane 
solution 


^0 HM 


number of modes sought o 


5o AMYK 


code calling for the use of a lumped or 
distributed mass model in analysiSo 


6o 3D 


code calling for inclusion or neglect of 
shear deflection ef f ects o 


7o HI 


code calling for Inclusion or neglect of 
rotational inertia effects » 


8 0 OMEGAG 


upper limit of frequency range of investi- 
gation (radlans/seco ) 0 


9 0 OMEGAI 


lower limit of frequency range of investi- 
gation and initial interatlon frequency « 

If a zero is placed in this field the 
search for mode frequencies will begin at 
oOl radians per second o 


10 0 GDaBC 


code calling for print out of the frequency 
and value of the frequency determinant for 
each iteration, the number of subsections 
assigned each component treated by lumping 
mass, the boundary condition state vector 
for both ends of the main member and of each 
branch, the number of components in the sys- 
tem, the frequency range of Investigation, 
the highest mode sought, and all of the input 
datao 
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Codes used in the first data card are as follows 



a. I0P=0 


both in- and out-of-plane mode fre- 
quencies sought. 


b. I0P=1 


In-plane frequencies sought. 


c. I0P=2 


out-of-plane mode frequencies sought. 


d . AMIK=1 


lumped mass model desired. 


e . AWYK=0 


distributed mass model desired. 


0 

II 

Q 


consider shear deflection. 


g. SD=1 


neglect shear deflection. 


0 

II 

M 


consider rotational Inertia. 


i. RI =1 


neglect rotational Inertia. 


j . GDA3C=0 


print solution only. 


k. GDaBC=l 


print solution and input data. 


1 . GDa3C=2 


print solution, input data, and 
graph data. 



B .72 Boundary Conditions 

Boundary conditions for the in- and out-of-plane case 
are read in separately. Codes outlined in B. 6 I apply to 
both cases. When both in- and out-of-plane mode frequencies 
are sought, the boundary conditions for the in-plane cases 
are read in first. 

The field specification for the first boundary condi- 
tion data card is: 

22F 3.0 

The first N+2 fields of this card must be filled where 
N is the number of branches in the system. The coded repre- 
sentation of the main member's left and right end boundary 
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conditions are entered In the first and last fields used 
respectively o The codification of the far end boundary 
condition for each branch is entered in the intermediate 
fields in accordance with the branch sequence number (Bo53)® 
When a boundary condition is not descrlbable by one of the 
codes given in Bo 6 l, a Oo must be entered In its data field. 
For each zero entered in the first boundary condition data 
card an additional data card is required on vihich is entered 
the actual state vector as given In Bo 62 or B 063 , codified 
in accordance with the last paragraph of B 06 I 0 These cards^ 
are placed immediately behind the first boundary condition 
data card their order dictated by the sequence number of 
the branch to which they apply (Bo53)o 
B .73 Branch Joinlnpc Anprle 

When the system includes branches, the angles at which 
the branches Join the main member are entered in one or two 
data cards as required whose field specification iss 
20 F 7o3 

For an N branch system the first N fields are used. 

The angles, in degrees, are entered in successive fields 
in accordance with the "sequence number (Bo53) of the branch 
to which they apply. 

B . 7^ Component Extensive Properties and Branch Component 

Identification 

A data card is required for each system component. 

The order of these cards is dictated by their component’s 

^Field specification 6F2 o0o 
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sequence nua.ber (Bo5l)= i'he card field specification iss 
3F 5o2 , 2F 6 c2 , 6F ?„2 , 12 

Field names and their data content are given, in order, 
below? 

D outside diameter ( inches; 

2 

DI wall thickness (inches), 

RHO radius of curvature with sign 

(inches) c Notes Applies to curved com- 
ponents only, A zero (0.) is entered in 
all other cases. 



TL length (Inches). Notes applies to 

straight components only, A zero (0,) is 
entered for all other cases. 



THBTA are central angle (degrees). Notes Curved 
components only, A zero (Ct) is entered in 
all other cases. 



CLX 
CLf 
CLZ L 
CTX 
CTY 
CTZ J 



spring constants 



(B.45^ 



NImBR branch component identification number (B,42), 
B,75 Intensive Properties of Components 

The first intensive properties data card has four 
fields whose specifications are 

F 7o3 , 2E 10,2 , F 2,0 



1-2 

These fields are left blank in cards pertaining to 
spring hangers, 

3 These fields are left blank except in cards per- 
taining to spring hangers. 
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The field names and their data content are given below= 

AAMU density of first component (IbSo/cu, ft,). 

AE elastic modulus of first component (psi), 

AG shear modulus of first component (psi), 

DD code, see ^ below, 

^ When the intensive properties of all components (spring 
hangers excluded) are identical, the DD field is left blank 
and the first intensive properties data card becomes the 
only intensive properties data card. If DD is non-zero, 
the data on the first card is interpreted as pertaining to 
the first component only and an intensive properties data 
card must be prepared for each of the remaining components, 

i.e,, the second through the last, spring hangers included. 
These additional cards are Identical to the first card in 
format and field content except that the DD field is deleted, 
A blank card must be inserted for each hanger. Intensive 
properties data cards are ordered in accordance with their 
component’s sequence number (B,51)<> 

B,8 Data Deck Assembly 

The data deck is ordered as follows: 

1, problem statement card, 

2, boundary condition data card(s), 

3, branch Joining angle data card(s), (deleted when 
system has no branches), 

4, component extensive properties data card(s)o 

5« component Intensive properties data card(s). 
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I 




The user may obtain solutions for any number of prob- 
lems desired by "stacking" assembled problem data card 
decks at the end of the program deck. Note that a blank 
card must be inserted between the last card of the program 
deck and the problem statement card of the first problemo 
A blank card must also be placed behind the last intensive 
properties data card of the last problem. DO NOT Insert 
blank cards between problems. The problem statement card 
of the N+first problem must be immediately behind the last 
intensive properties data card of the problem. 

3.9 Flexibilities of the Program VIPIPE 
B.9I Shear Area Form Factor 

A shear area form factor of one half is built into 
the program, [s] It may be altered by changing card 000580. 
3,92 End Conditions 

A properly constituted boundary condition state vector 
must have three zero and three non-zero elements. Fre- 
quently, however, an actual end condition may not meet 
this specification, e.g., an end of the main member or of 
a branch may be supported, or connected to a piece of 
machinery in such a manner as to approximate being held 
in a rubber mounting block (Fig, 3.92-1). To program this 




Fi', 92-I Fictitious end condition 
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end condition, the user should make a good approximation 
of the elastic properties of the mounting, support, or 

connection in terms of the linear and torsional spring 

2 

rates of the program (B,45, then attach a fictitious, 

short (perhaps *01 in*) massless straight section to the 
end under consideration, and specify in the boundary con- 
dition for that end that the end is free* 

Through the use of fictions of this sort, properly 
constituted end condition state vectors may be constructed 
for any real situation v/lth little or no loss of accuracy* 

B* 93 Subsectioninp; 

Subsectioning for the lumped mass treatment of com- 
ponents is carried on in subroutine SUBSEC^* The subroutine 
is set up to divide a component into a maximum of ten sub- 
sections, The writer, through the analyses of several 
systems, found that this figure gave good accuracy of 
solution for a reasonable expenditure of machine time* 

For all components having a length-to-diameter ratio 
greater than six, the number of subsections is computed 
as two times the value of a parameter HM which appears 
in the subroutine. The user may control its value by 
altering card 004?20 in the subroutine* 

^See App, A, sections A*22, and A*32* 

^i,e,, treat the mount as a hanger* 
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B.9^ Values of the Frequency Determinant and Solution 
iicceptability Criterion 

The frequencies for V'Jhlch the value of the frequency 
determinant is zero are the natural frequencies of the 
system which it describes » In light of the number of sig- 
nificant figures available in a computer and of the effects 
of round-off, however, insistence that iteration and con- 
vergence operations continue until the value of the 
frequency determinant actually became zero would require 
an excessive expenditure of machine time* Instead, itera- 
tion proceeds until the sign of the value of the frequency 
determinant changes, then a convergence routine is carried 
out until the difference betv^een the values of frequency, 
for which the corresponding value of the frequency determinant 
have opposite signs, is less than a specified value^CR* 

The result of a linear interpolation between those fre- 
quencies is accepted as a natural frequency of the system. 

The specified frequency difference CH is termed the solution 
acceptability criterion. 

To compute a starting frequency increment, the program 
constructs a synthetic straight pipe, equal in length to 
the length of the main member, and having diameter, wall 
thickness, and intensive properties equal in magnitude to 
a weighted average of those properties of the composite main 
member. The fundamental frequency of the synthetic pipe 
is computed, and one half of its value is taken as the 
starting frequency increment!” The initial solution 

^This increment may be changed by altering card OOIO 9 O. 
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acceptability criterion is taken as one-twenty-fifth of 
the starting frequency increment o The initial solution 
acceptability criterion may be changed by altering card 001110» 
The starting frequency increment and the initial solu- 
tion acceptability criterion are used until the second 
natural frequency is found, after which the frequency 
increment is recomputed, as it is after each succeeding 
natural frequency is determined, and given a value equal 
to one half of the difference between the last two fre- 
quencies found. The new solution acceptability criterion 
is then recomputed as one-fiftieth of each new frequency 
increment. The new solution acceptability criterion may 
be changed by altering card 002800, 

B,95 Flanges and Couplings 

Pipe couplings and flanges represent real "lumps” of 
mass mounted in a piping system. They may be treated as 
individual system components. If this is desired, a com- 
ponent sequence number and a branch component identification 
number should be assigned in the usual manner (B,51, and 
B,52), and appropriate extensive and intensive properties 
data cards (B,?^, and B,75) should be included in the data 
deck. These components are treated in the program, as 
rigid sections having distributed mass, 

B»96 Hangers Exhibiting Coupling 

Hangers such as that shown in Fig, B,96-1 can be 
treated vjith the tools available in VIPIPE by considering 
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the mounting arm as a branch having a hanger 




Fig, B.96-I Cantilevered hanger 
and a free end. Recourse to similar methods will permit 
treatment of almost any hanger configuration. 

B .97 Systems Having Unusual Flexibility 

When the main member of a system is essentially 
cantilevered the computed starting frequency Increment may 
be so large that the first two system natural frequencies 
are not found. For analysis of such systems It Is recom- 
mended that the constant 4.7300^019 in card 001080 be 
changed to I. 875 IOO 67 . Similarly, If It is suspected that 
the fundamental frequency of the system is less than 
.01 rad. /sec. a suitably lower starting frequency should 
be programed (B. 7 I). 
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APPENDIX C 



FLOW DIAGRAMS 



C«1 General Remarks 

Flow diagrams are included for all parts of the main 
body, and for several subroutines. Those subroutines for 
which no flow charts are Included, with the exception of 
INVERT, Incorporate, at most, two decisions in their 
structure. As applicable they sense on SD and/or RI ( A.31) 
to change the value of one or two elements of the transfer 
matrix constructed by the particular subroutine. 

Subroutine INVERT performs a matrix Inversion by a 
pivotal technique. Sufficient comment statements are 
Included in the subroutine to define its structure. If the 
user’s FORTRAN library includes a satisfactory single 
precision matrix Inversion subroutine, INVERT may be deleted 
after suitably modifying the dimension (cards OIOO 7 O, 

010080 BRANCH, and 011080, OIIO 9 O BRANCO) and call (card 
010230 BRANCH, and OII 66 O BRANCO) statements of subroutines 
branch and BRANCO, 

C.2 Index of Flow Diagrams 
C.21 Main Body 

Flow Diagram Page 



INPUT 

INVARIANTS 

CONTROL IN-PLANE i 
CONTROL OUT-OF-PLANE 
ITERATION 
OUTPUT 



C-6 

C-6 

C-7 

C-9 




C-1 



m 




C.22 Subroutines 



Flow Diagram 



Page 



SUBSFC 

IvIilTMUL 

BRANCH 

BRANCO 

STaVEC 

STAVEO 

FINnAT 

FINBRa 

DELilA 



C-10 
C- 8 
C-11 
C-12 

C -13 

C-13 

C-14 

C-14 

C-14 



C.3 Flow Dlafrrams^ 



C-3 through c-l4 



Note: Since CONTROL IN-PLANE and CONTROL OUT-OF-PLANE 

are Identical In structure, only the former flow diagram Is 
Included. They differ only in the subroutines called and the 
state vectors used (A.22). Subroutines SUBSEC, MATMUL, INVERT, 
FINMAT, FINBRa, and DELMA are common to both cases. 



^A complete list of symbols is given In Appendix A, 
section A. 3 . 
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JJPTTT hJjO W PXA.^I- AM ) 






AMYK = AWYK-1. 

OME^AO =OME''AI 
FEAC= .{> 

L = 1 
M =1 

AA- PP = CC = r,0 = EE=0. 
AaMH = AL -AD=AID=0. 
KK -0 

NVN = NPP +-1 



lHETR(l) = THS'rA(l) 

DT(I) = Dif I) 

RAMir( I ) = AOT ( I ) 

DI(I) = DI(I) - Z.O#m(I) 
1=1, NS 



I (;0 TO INVA RIA NTS 



INVARIANTS FLOW DIA >HAM 




C-'l 



TM' AT.iAW 'i’a ■'•'j.O’v 




C-5 



■■"'I, tn Pr,Af;E Flow FTAilPAM 




c-c. 




1 



ITHPATION 



FT,OW DIA 'hMA 




STBROnTINE KATMHL FI.OW DI-ti'rRAM 





r-S 



pj; TPy f y 




|.TJ= 



IP(I) 


= 0,1 


Z(I) 


II 

O 


II 


,MS 1 







[input. 



>0 









MN = AA=BB=CC =GG=-KK=0. 
ADOM - BDOM 
CR = BR 

OMEGA ( 1 ) = OMEGA I 
OMEGA 0 = OMEGA I 
DOM = ADOM 






IGO TO Addrep 

[control out- 


s 998 ofl 

of-plane| 



C-3 



OlLTi - . BJ’PSKC DIA . rJVW - 



■'v ’ OD=APS.'-'(- ofl 




C-ID 



■'kU 



ctADCH. PLOv; CIj, 



T 

po = l) 




r - ’ -1 

1.1 = ] ^ ’<=i,6j 

. j;: ... 

I’r( 4 ^ 4 )= 'r(6.b)= TT{^),6)= 1. 

^T( 4 ,’-,)= RH ,T.) 
TI(5,L)= R(S,L) 

(^(6,r,) = H(3,L) 

• L L=J .3 



rTKr,.n^ 



I R(J,K)= 0 . I 
|j= 1,5 K= 1 . 5 | 



|R(J,K)= R(J,K)+R 2 (J,L)*R 1 (L,K) 

I L= 1. 5 K=l ,3 J = 1 .3 I 



I SP= SINF(PHI) 

I CP= f-O SF(P Hl)l 



oordinate | 
01 ar-d r, 2 j 



I R1TJ,K)= 0. I 
IT = 1, 5 K- 1,5 I 



I R1(.T,K)= Rl{.T,K)+d2(.T,T,)^^R(L,K) I 

_L=1,3 K=l,3 .T=l,3 

^ 

I T = 1,3 K = 1.31 
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S:'pRqUTI!'E._PPA!'._CO KLOV/ DIA. IRAM 



. i . 



Peordar bour.dbry conditlop state 
vector such that the first three 
elemet'ts pertain to displacements, 
_the, last three to. forces 



n^eorder branch transfer matrlxl 

L to make it compatible with the 
reordered state vector,. _ | 








1 R1(1-3,K)= 
R2(1-3,K) = 
1 K=1 


UV(1-3,K) 

UV(4-6,K) 

3 







1 CALL INVERT I 
I (Invert R1 ) I 



n-'i' - , 

L'I= 1^3 K= 1,3| 



S ( J , K ) = ■ S ( ■/, K )+ R 1 ( j , L ) * G1 ( L , kTI 

- ■ ‘ 3_J 



L= ..K = 1,3 



U(l,l)= U(2,2)=U(3,3) = 1. - 

U(4,4)= U(5,5) = U(6,6) = 1. 

U(2,l)= S(l,l) U(b,4) = S(2,3) 

U(2,5)= S(l,2) U(6,l) = S(3,l) 

U(2,4)= S(l,3) U{6,3)= S(3,2) 

U(5,l)= S(2,l) U(6,4)= S(3y3) 

U(5,3)= S{2,2) 



r~ etorn1 



rR( J,K)= R( J,KH-R2( J,L)*R1(L,K) 

1 L=l,3 K = l,5 J=l,3 I 



|SP= SINF(PHI)|' 
|CP = C 0 SF(PHI)| 



, Il_ 

I R1(J,K)=0. 
|J = 1.5 K =l,: 



K = li3 



C-l^ 



__S^TCPO'rTTWK S'l^AXTC PLOW . PIA . 




( 4 ,N ) = EC ( fa,N ') =* BC ( 6 ,n') = 1, I — ^ »- f RETU RN | 



^}Cg,N )- RC(2,N y=Br (3',N T= A >- {^'nJR^ 



• ] PC( B.~N 1 = BC (5 tN ) r . P? ,W )=iyi >- j RE^'-'pil 



■j' PCd .N ) = P<'(3,P ; =iBC(fi.N ) = l.| =Hb BTURn1 

( ®i^TLP^j'?JLN.r=' PC ( 6 ",nTVi 



SUBROUTINE STAVEO FLOVJ DTA'-BAM 



Hpc(2 . n)^c(5 .n)=rc(6 .nT^iTI H BEriRN I 

■rBC(l,N)= BC(3.H)= B C(4,N)« l. | >- | RETURN 1 

PC ( 8 ,N H~BC~( 47i^ BcTMT ) = 1 . J H PErTRN~l 



- - ■■— > |rv.t tikh I 
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^UPI^PJITINE . F_3;;KAT tLOVii DIA-;hAM_ 



i-^ " ^ V ( .f , k') = u ( .f , k ) 1 

■ L K= 1,6 J 



Vf.J,K)- 0. r- 

^.I -3-1,6 K=-l,y 



H.^ETURNH 



V(J,K)- V(J,K)+u'(J,Lj»UU(L,KT]-. -^RS TOhjil 

.L=l,6 K = 1,6 J- •J 



<wr 



v(j,kj*u(j,k) 

■T = 1, 6 K = l,6 



■ j'p.lTU Rk'i 



I V( J,K)= 0. = V( J,K>U( J,L)*VV(L,K) F 

I J - 1.6 K =1,6 I , I L=l ,6 K = 1.6 J= 1.6 J 



S^n«HOnTIHE ..mWA- .ZLQ'N DIAGRAM 
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APPENDIX D 



PROGRAM LISTING 



D.l Index 



Dell Main Body 



Section 


Sequence Numbers 


Page 


INPUT 


000040-000760 


D-2 


INVARIANTS 


000790-001170 


D-3 


CONTROL IN-PLANE 


001200-001810 


D-4 


CONTROL OUT-OF-PLANE 


001840-002450 


D-5 


ITERATION 


002470-003260 


D-6 


OUTPUT 


003280-004660 


D-8 


SUBROUTINES 


SUBSEC 


004690-005140 


D-10 


DISTM 


005170-005710 


D-11 


DISTMO 


005740-006280 


D-12 


SFIELD 


006310-006640 


D-13 


SFIELO 


006670-006910 


D-14 


C FI ELD 


006940-007370 


D-14 


C FIELD 


007400-007890 


D-15 


POINT 


OO792O-OO8O9O 


D-16 


POINO 


OO8I2O-OO83IO 


D-17 


MATMUL 


008340-008900 


D-17 


RIGID 


008670-008640 


D-18 


RIDIO 


008930-009160 


D-18 


STIFCO 


009190-009410 


D-18 


STIFOO 


009440-009660 


D-19 


HANGER 


009700-009850 


D-19 


HANGED 


009880-010030 


D-20 


BRANCH 


010060-010730 


D-20 


STAVEC 


010760-011040 


D-21 


BRANCO 


011070-012210 


D-22 


STAVEO 


012240-012480 


D-24 


INVERT 


012510-013220 


D-25 


FINMAT 


013250-013380 


D-26 


FINBRA 


013410-013540 


D-2 6 


DELMA 


013570-013750 


D-27 



D.2 L^stl^ D„2 through D-2? 



PRCCR/'w VIPIPl 



c 

C INPUl 

C 

c 

NP = 0 
JJ=1 

M500DIMLNSICN aj((5C>,AJ(50),R(t;C),AIX(S0),AIY(50),PHI12C),SVI(t),1), 
1SVB(6),1;(50),CI(5C),TL(5C),/'. (50) , G ( 5 0 ) , E ( 50 ) , AML' ( 5C ) , THE TA ( SO ) , RHO 
2(50),CMLCAM(20) ,CEL(5CO),rVtCA(5CC), A(6,6(,P(A,<.),U(t,61,CLX(5C),C 
3L¥(5C) ,CLZ(5C),CTX(501,CTY(S0),C1Z(5C),V(6,6),VV(6,6) ,LU(6,6),SVL( 
i46,l),SV(22),SVIP(A,22),SVP(Z2),SVnP(6,2> l,NNBR( 50) ,BAMU(-SC I ,THETR( 
520) ,DT(5C),P(20) 

REAL 15CC,NS,NRR, IGP.HM, AMVK ,SD,RI ,TMECAG,OKEGAI,GDAEC 

1500 FCRyAT(3I2,l(E3.C,2Fe.O,F2.C) 

NyBR=NBR+2 

IF( ICP-1 ) 1501 , 1501, 1506 

1501 REAC 1502, (svm, 1 = 1, ^yeR) 

1502 FrRMAT(22F3.0 ) 
cc 1505 j=i,NyeR 
IF(SV(J) 11503, 1'03, 1505 

1503 READ 150»i,(SVIP( I,J), 1 = 1,6) 

1504 FCRMAT(tF2.0) 

1505 CCNTlNUf 

IF(IOP) 15C6, 1506, 1510 

1506 REAL) 1507, (SVC( 1 )', 1 = 1, NMER) 

1507 rCRMAT(12F3.0 ) 

CC 1510 J=l,NyBR 
IF(SVC(J) 11508, 1500,1510 

1508 REAL 1509 ( S VCP ( I , J ) , I = 1 , 6 ) 

1509 FrRMAT(6F2.0) 

1510 CONTINUE 

IF(NeR)15 13,15H,1511 

1511 RFAD 1512, (PHK t ), 1 = 1, N6R) 

1512 FORMAT ( 1CF7.3) 

15130REA0 1514,( C( I ) ,CI ( I),RHC(I ),TL( I ) ,THETA(I ) ,CLX(T ) ,CLY(I ) ,CLZ(I ) ,C 
1 TX( I ),CTY(I ) ,CTZ ( I 1 ,NNPR (I ) , 1 = 1 ,NS) 

1514 FCRMAT(iF5.2,2F6.2,6F7.2, 12) 

READ 1515 , AAMU, AE, AG,CC 

1515 FORMAT (F7.3,2E10.2,F2.C) 

IF(CC) 1516,1516, 1518 

1516 or 1517 1=1, NS 
AMU( I )=AAMU 
E(I)=AE 

1517 G( I)=AG 

GO 1C 1520 

1518 AMU( 1 )=AAMU 
EM )=AE 
G(1)=AG 

READ 1519, ( AMU( I ) ,E( I ) ,G( I ) , 1 = 2,NS) 

1519 FORMAT ( F? . 3 , 2E 1 0 . 2 ) 

IFCCMEGAI ) 1520, 1520,1521 

1520 CMEGA1=.01 

1521 CMEGAM )=CMtGAI 



0C0C2O 
0CC030 
00CC4C 
0CCC5C 
OCOC60 
0CCC70 
OOOCRC 
OC0C90 
.000100 
oconc 
000120 
OC013C 
000140 
0001 50 
000I6C 
000170 
000180 
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AyYK=ANYK-l .0 
CMCGAC = C^•EGAI 
FF At = .‘^ 

L = 1 
N=1 

AA=0.0 

BP=0.0 

CC=0.0 

GG=0.0 

KK=C 

AAMC=0.0 

EE=0.0 

AL=0.0 



AID=C.O 
NNN = NBR+ 1 
DG 1522 1=1, NS 
DT« I 1=111 ( I) 

BANU( I l = ArU( I I 
THETR1I)=THETAI I) 

1522 DI ( I l=C( 1 )-2.C*Cl( I ) 

C 

C 

C INVARIANTS 

C 

DO 1C8 1=1, NS 

IF(CLX(Il+CLY(II+CLZ(I)*CTXm+CTYlI)+CTZ(Il)10e,101,lC8 

101 AMUI I )=AVU( I ) *3. 1415926 5i^*(t ( 1 1 * •2-lH ( I ) * • 2 ) / ( 69 1 2 . * 32. 1 7* 1 2 . ) 
IF(RHC(iniC2,103,102 

102 THElAm=THLTAn)»3.1U15S2654Zl80. 

TL( I )=THETA( I )*ABSF(RHC( I ) 1 

103 AJT( I 1=3. 141 59265U*(D( II • »4-C I ( I ) * »4 ) / 32 . 

AJ( I 1=AJT( I 1/2. 

R( 1 )=1 .0/(AJ( n»L( I) ) 

AIX( 1 )=SGRTF( (C( 1)«*2 + DI ( I 1«*2)/F. ) 

AIY( I l=AIX( I ) /I .414214 
IF(AVYK) 104, 1C5, 1C4 

104 Z( I ) = 1 .0 
IF(RHOm)l05,106,105 

105 HHM=HV 

CALL SLISSEC I RHC ( I ) , THE T A ( II , D ( I ) , TL ( I 1 , , Z ( I ) 1 

106 IFINNPRI I ) ) 108, 107, 108 

107 AL = AL + TL(n 
AO=AC + D( n*TL I I ) 

EE=EE+EI I )*TL ( I ) 

. AID=AIC + DI(II*TLm 
aahl;=aamu+ai-lm I ) *TL( I l 

108 CONTINUE 
AD=AC/AL 
AIO=AIC/AL 

AFR = 3. 14 159 26 54Z64.*(AC»«4-Ain**4) 

FR=4.730C4019«*2«SCRTF(EE*AFR/(AAMU*AL»»4J1 

AU0P=FR/2.0 
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0C073C 
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000760 

000770 
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000810 
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000830 

000840 

000850 

000860 

000870 

000880 

000890 
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000910 
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000950 
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l)Cy = AL'L;- 
C'»=FP/15.C 
BDON = ACOF' 
eR = LR 

cr: ICO 11 = 1 , NTR 

109 PHI ( 1 I )=PHI ( I I 1 »3. 1M 1 5926 54/ 180.0 
IFt ICP-1 ) 1 1C, 990, 996 

110 MN=1 
C 

C 

C CONTROL IN PLANE 

C 

999 IF( SV( 1 1 I 101 1, 1C21 , 1C1 1 
1C11 CALL STAVLC I S V ( 1 1 , 1 , 5 V I I 
GO TO 1041 
1021 or 1C 31 J=1,6 
1031 SVI ( J, 1 )=SVIP( J, 1 1 
ICUl IFISVINFPR) ) 1051 , 1061, 1C51 
1051 CALL STAVECISVIN^PR), 1 ,SVC) 

GO TL 1CB1 
1061 00 1071 J=l,6 
1071 SVE( J, 1 1=SVIP(J,NPPR1 
1081 JK=C 

IF (NCR 11000, 1C0C, 108 5 
1035 or 1087 N=2,NNN , 

IFISVINl 1 1087,1087, 10R6 

1086 CALL STAVECISVIM ,N,SVIP 1 

1087 CONI INUfc 
1000 N=2 

11 = 1 

00 1001 1=1, NS 

P( I )=CLX( II +CLYI 1 I+CT2 ( 1 1 

IFIPin 137,36,37 

37 CALL HANGERICLXI I 1,CLV( I I,C1ZI II ,UI 
GO TC 62 

38 IFIAPYK139,41,39 

39 IFIRHOI I 1 140,42,40 

40 IFIZ I I 1 143, 45,45 

41 IFIRHCI 1 1140,48,40 

42 1FIZII1146,47,46 

43 CALL CFIELD I R I 1 1 , Z ( I 1 , G ( I I , SO , RPC ( 1 1 , TH ET A ( I 1 , 0 ( I ) , D 1 1 1 1 , E (II , A ) 

44 CALL POINT ( ARU ( 1 1 , A 1 Y ( 1 1 , O' EG A ( L 1 ,T L ( 1 1 , Z ( 1 1 , R I , R 1 
CALL HATMUL ( A,P,Z( I.),U1 

GO rc 62 

45 CALL STIFCO I THETA ( I ) , RHC ( I 1 ,U 1 
GO TO 62 

4 60CALL DISTF* ( ANU ( 1 1 , A I Y (II ,0 P EG A ( L 1 ,T L ( 1 1 ,R ( I I , 0 ( 1 1 , D 1 1 1 1 , G ( I I , SC , R 
1 I ,E( I 1 ,FFAC,U I 
GO rc 62 

47 CALL RIGID I AKU ( I I , TL ( I 1 ,CMcGA ( L ) , AI Y ( I 1 ,RI , LI 
GO TC 62 

43 IFI Z( 1 1 149,47,49 

49 CALL SFIELD ( C I ( I 1 , TL ( I 1 , R ( I 1 , Z ( 1 1 ,3 ( I I , SO , 0 ( I) , E ( I I , FFA C , A 1 
GO rc 44 
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001420 
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001450 
001460 
001470 
001460 
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001500 
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001520 
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001560 
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001580 
001590 
001600 
0016 10 
001620 
001630 
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62 1 F (NNPR( I ) )65,6i,65 

63 CALL FINPAT (L,LL,I,V) 

DC ok J=l,6 

DC 64 K=1 ,6 
6U L'Lit J,K)=V(J,K) 

GO IC 10C1, 

65 CALL Fir,?RA ( L, VV, NNCR ( I ) , V ) 

CO 66 J=1 ,6 

DC 66 K=U6 

66 VV( J,K)=V(J.K) 

IF(^^SR( I) ) 68,67,67 

67 lF(NNbR( 1 )-3ll001 ,68,68 

68 CALL BRANCH(SVIP,VV,PHI( II),N,U) 

N=N+1 

n = ii + i 
GO TC 65 

1001 CONTINUt 
GO TC 2000 

C 

C 

C CONTROL CUT OF PLANt 

C 

998 IF(SVC(1 ) 11013, 1023,1013 
1013 CALL STAVEO(SVCm,l,SVI ) 

GO TO 1043 ' ' 

1023 DO 1033 J=l,6 
1033 SVK J,1 )=SVCP(J, 1 1 
1043 IF(bVO(NKBR)) 1053, 1C63, 1C53 
1053 CALL STAVE0(SV0(NFBR1, 1,SVE) 

GO TC 1083 
1063 DO 1073 J=l,6 
1073 SVE( J, 1 )=SVCP( J,NPBR1 
1083 JK=1 

IF(NF.Rn0C2,lC02, 1093 

1093 DO 1095 N=2,nKn 
IF(SVO(N))1095, 1095,1094 

1094 CALL STAVEO(SVC(N),N,SVCF) 

1095 CONTINUE 

1002 N=2 
11 = 1 

DO 1003 1=1, NS 
p( n=cTxm+CLZ( ii+cty( i i 
1F(P(1 111,2,1 

1 CALL HANGEO(CTX( n.CLZl I 1,CTY( II ,U) 

GO TC 14 

2 IF(APYKl3,n,3 

3 1F(RHC(1114,8,4 

4 IF(Z( I 116,5,6 

5 CALL STIFCO ( THE 1 A (I | , RHC ( I 1 , U 1 
GO 1C 14 

60CALL CFIELOI AJT( I) ,R(11, Z( I 1,G( I 1 , SD,RHG(II , THETA( I 1 ,0(1 1 ,01 ( I 1 ,FF 
1 AC,A1 

7 CALL PCINC: ( APU( I 1 , AIX( I 1, AIY( I 1 ,CMEGA(L 1,TL( I 1 ,Z( I 1 ,RI, 81 
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cr Tr \u 

8 IF(Z( I) ) 10,9,10 

9 CALL RICIC ( ANU( I ) , IL ( 1 ) , AIx( 1 ) , ( VEGA (L) , AI Y( I ) , AJT( I ) ,R I ,U) 

GO rc 1 U 

lOOCALL OISTM](AyU(I),AIX(I),AlY(Il,OMEGA(L>,TL(n,AJT(I),R(I),L(I),C 
1 I ( n ,G( I ) ,SD,RI ,FFAC,li) 

GF rc lA 

n IF(RHC( I) )U, 12,u 

12 I F ( 2 ( 1 ) ) 1 3 , 9 , 1 3 

13 CALL SFICLO ( C I ( I 1 , TL ( I ) , A J T ( I ) , R ( I) , Z ( 1 ) , G ( I 1 , SO , C ( I ) ,F F AC , A ) 

GO TC 7 

14 IF(NKhR( in 17, IS, 17 

15 CALL FINMAT (U,GL, I, VI 
DO 16 J=l,6 

DO 16 K=l,6 

16 UU(J,K)=V(J,K) 

CO IC 1003 

17 CALL FINPRA ( L , VV , NNER ( I ) , V ) 

DO 18 J=l,6 

DC 18 K=l,6 

18 VV( J,K)=VU,K 1 
IF(NNBR( 1)182,19,19 

19 IF(NNFR( n-3) 1CC3,02,R? 

82 CALL BRANCO(SVCP,VV,Ph'I( II),N,IJ) 

N=N+1 
11*11+1 
GC 10 15 
1003 CCNTINUL 
C 

C ITlRATION 

c 

2000 CALL DELPA(UU,SVI,SVE,V) 

ODEL(L)=V( 1, 1 )*V(2,2)*V(3,3)-VM , 3)*V ( 2,2)»V(3, 1 )+V( 1 ,2)» V(2,3)»V( 3 
1,1)-V(1,2)*V(2,1)*V(3,3)+V(1,^)*V(3, 2)*V(2,n-V(l,ll»V(3,2)*V(2,3) 
GO TO 61 

21 PDEL=OEL(L) 

CO TC 25 

22 IF( AA)23,23,24 

23 PDEL = DEL(L-1 1 
GO TO 25 

24 PDEL=CELC 

25 IF(PCEL*0EL(L))26,26,27 

26 IF(AA)28,28,57 

27 BB=0 

IF( AA135,35,52 

28 AA=1.0 
. BB=1.0 

DELC=LEL(L-1 ) 

CPEGAC=CyEGA( L-1 ) 

GC TC 59 

2900MECAy(y)=aFEGAC+(0yEGA( D-CMEGAC) «A 8 SF ( CELC ) / ( ABSF ( DELC ) ♦ 

1 A8SF(CEL(L) ) ) 

GC TC 32 

30 CMEGA(L+1 )=0MEGA(L )-.9*ury 
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31 L=L+1 

IF( JKl lOOC, 1000, 1002 

32 BM=M 

ir(HF-En3DOO,3CCC, a 

33 y=w+i 

IF(N-3)3U0,5i4,3U ' 

34 ADCM= (QN'EGAN'CN-l 1 -0^ EG A'* ( r-2 ) ) • . 

CR = ACC»'/bC. 

340 AA=0.0 
BP = 0 

cc=c 

CG=C.O 
DCM = ADOf' 

DELC=0 

O.MECA(L+l)=CMEGAy(M-l)+2.»CR 
CMEGAC=Of'EGA( L+ 1 1 
GO TC 31 

35 1F(lKECA(L)-C^EGAG)55,60,6C 

50 0^'EGA(L•^1 ) = CMEGA( L»-(0VEGA(L l-CNfCAC 1 /?.C 
IF(Lr'EGA(L)-0^£GA(L + l 1 -CR ) 5 1 , 5 1 , 3 1 

51 CC=1.C 
GO TC 51 

52 DELC = UEUL) 

CF'EGAC=OMEGA( L) , 

53 CMEGA(L+1 )=CNEGA(L)+D0N 
GO TO 51 

54 IF(GELC*0EL(Ln29,2q,55 

550CMEGAHM)=0MbGA(U + (CMEGA(L-l l-OfEGA (L) 1*ABSF(DEL(LI ) / ( A BSF ( OEL ( L 1 
1 ) + ABSF(DEL( L-1 ) ) ) 

GO TO 32 

56 IFICKEGA(L)-0I-EGAG)21,21 ,22 

57 IF(BB)5B,5a,5C 

58 Pfi=1.0 

59 IF ( OMEGA! L)-QPEGA( L-1 ) -CR » 29 , 29 , 30 

60 Fi=M-l 

60 TO 3000 

61 IF(GG)70,70,72 

70 IFIDEUL) 177,71 ,77 

71 GG=1.0 

0MEGA!L+1 l = CMEGA(U-CR 
GO TC 31 

72 IF(GG-1.0)73,73,7U . 

73 GG=2.0 
HOEL=CEL(L) 

CMEGA(L+n=CMEGA(Ll+?.*CR 
. GO TO 31 

74 TDEL=DEL(L)»HCEL 
IF(TCEL)75,76,76 

75 Cf'EGAM(MI=CMEGA(L-2) 

GO TC 33 

76 KK=1 
M = M-1 

GO TC 30CC 



0C272C 

002730 

0C2740 

0C2750 

0027C:C 

002770 

0C27PC 

0C2790 

Q02ECC 

002810 

002820 

0C2Q3C 

002840 

OC2B50 

0C2560 

002670 

002SrC 

00289C 

0C2900 

002910 

002920 

002930 

0C2940 

002950 

00296C 

002970 

002980 

0C299C 

OC3COO 

003010 

OC302O 

OC303C 

OC3C40 

003050 

003060 

003C7C 

0C3C80 

003090 

003100 

003110 

0C3120 

003130 

003140 

003150 

003160 

003170 

0C3180 

003190 

0032CO 

003210 

003220 

003230 

003240 

003250 



D-7 



77 IF(CC)56,56, 'U 00526C 

C 0C3270 

C OUTPUT OC32GC 

C 0032'5C 

3000 IF(NP1UOO,UOO,UC! 0C33CC 

400 PRIM 401 0C3310 

401 FORMAT (IHI//) 0C3a2C 

PRU.T 4C2 003330 

4020FCR.''AT {15H PRCOR am' v I P I PL , - 5X , 1 7HG. E. FINK NH A- 1 , 5? X, 1 3HF cPRU AR 0C3340 
1Y 1964// , lOX, ICPH MnnAL FREQUENCIES OF A PLANAR PIPING 003350 

,1 SYSTLV ARE UET'RMNEC BY AN ITERATIVE PRCCFCURE USING T>E /29H PE 0C3 360 

3THCL CF TRANSFER N A TP I C E S . / / / M 9h ..a*..... 0033 7C 

003380 

5>. //) 0C3390 

403 IFINPI405E, 4050,404 0C340C 

404 PRINT 4C5 003410 

405 FCRPAT (IM//) 0C34E0 

4053 NM=i 0034i0 

4059 IF(JK)406C, 4060, 4065 OC3440 

4060 PRINT 4070 JJ 003U5C 

4C700FCRPAT I 5CX , 8HPRGPI £P, I 2 , // 3 5X , 4 7H IN PLANE NODE FREQUENCIES (RAOIA 003460 

INS PER SECOND) // 45 X , a hPCL E , 1 5 < , 1DH FREQUENCY /) 003470 

GC IC 4079 OO340C 

4065 PRINT 4075 JJ 003490 

4C750FriRPAT(5CX, BHPRrULEP I 2 , / / 3 X , 5 1 hCUT OF PLANE PCOE FREQUENCIES IRA 0C35C0 
IDIANS PER SLCCNC) / /4 5 X , 4HPCDE, 1 5 X , 1 OH FREQUENCY /) 003510 

4079 P.H = P 00352C 

PRINT 400 I N, OPEGAN. IP ) ,P= 1 , PH) 003530 

408 FCRP.AT ( 46X, I 2, 1 7X , F 1 4 .0 /) 003540 

IF(KK)4C80, 4090, 4080 OC355C 

4080 PRINT 4085 00356C 

40850FORPAT (10X,69HPRCPLEP TERPINATEC DUE TO SIGNIFICANT FIGURE LIPITAT 003570 

IION CF COMPUTER. /) 003500 

4090 IF(APYK)412,409,412 003590 

409 PRINT 411 003600 

UllOFORPAT (16X,53H A LUMPED MASS APPROACH WAS EMPLCYED, A 003610 

1N0 ) 0C3620 

GO TC 414 003630 

412 PRINT 413 003640 

413 FORMAT (24X,45HA DISTRIBUTEE PiASS APPROACH WAS EMPLOYED AND 1 003650 

414 IF(SD)415,U15,420 003660 

415 IFIRI >416,416,418 003670 

416 PRINT 417 003680 

4170F0RPAT (2X,67hTFE EFFECTS CF SHEAR JEFLECTICN AND ROTARY INERTIA W 0C3690 

lERE CONSICEREC. /I 003700 

GO TC 425 003710 

418 PRINT 419 003720 

4190F0RMAT (3X,63HTHt EFFECTS CF SHEAR UEFLECTICN WERE CONSIDERED WHIL 003730 

IE THOSE OF /3X, 35HR0TAT ICNAL INERTIA WERE NEGLECTED. /) 003740 

GO TC 425 0C375C 

420 1F(RI)421,421 ,423 003760 

421 PRINT 422 003770 

4220F0RPAT (3X,71HTHE EFFECTS OF ROTATIDNAL INERTIA WERE CONSIDERED WH 003780 

IILE THOSE OF SHEAR / 3X , 2 7HD E FLEC T I ON WERE NEGLECTED. /) 00379C 
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GC IC L 
423 PRINT 4^4 

ItnUOFGRFAT (;X,6rhTHF fFTECTS I'.F SHEA. i:LFLFtHn^' ANl: ROIARY INFP.TIA |«E 
IRE NECLtCTEl'. /) 

425 PRINT U4C 

4400 FCR RATI IIPH 

1 . * ) 

IF( GOAEC-1 . )4 5 1 ) , 4520, 44C 1 
4401 PRINT 4HC2 

44020FnRFAT( j 5H SECTICN LENGTH A.NL SU f> S EG T ION 0 A T A/ 2 i X , 1 4 L SEG T I GN NLRRE 

IP ,Px 5HLEN0TH rr sect iom inches ),:■<, 2 ihnurh e'» cf secthns ) 

PRINT 44 3 ( I,TL( I),2( I), 1=1 ,NS ) 

443 FORMAT ( 2 9X , I 3 , ^ 3 X , F7 . 2 , i 1 X , F4 . 0 » 

LR = L 

PRINT 441 

4410FCRRAT (11H GRAPH DATA / 2 "X , 1 6H I TERA T I CN NU^eER , 4X, 3CHFR CCUENCY IR 
lADIANS PEP SECCNC) , 3X, 35 h VALUE FRuGUCNCY DETERMINANT ) 

PRINT 444 (L,CRE0A(L1 ,DEL(L),L=1 ,LM) 

444 FCRRAT I 2 9X , I 3 , 2 3 X , F 1 1 . A , 1 G / , E 1 5 . ? ) 

4520 PRINT 442 

442 FORMAT ( 1 X, 10H INPUT DATA /I 
PRINT 4521 

45210fORHAT(1X,9HCCyPCNENT,4x,P,HCIAVETER, 5x, 14HWALL THICKNESS ,5X,6HLEN 
1GTH,5X,9HRACIUS OF , 9X , SH INCLUDED , 5 X, 7HCLNS I T Y , 5X , THE L AST I C , 5X , 5H SH 
2EAR /57x,9HCURVAtU'RL,5X, 12FANGLL CF ARC , 1 EX , 7HMCDULUS , 5X , 7HVCDULUS 

3 /) 

or 4525 1=1, NS 

IF(P( I 1 )4524,4524,452 3 

4524 PRINT 4522 I , C ( I ) , DT ( II , T L ( I ) , RHC ( II , THU TR ( I 1 , C AMU ( I ) , E( I) , G (I ) 

4 5220FDRl''AT (:3X,I3,Px,r7.:i,lOX,Fd.3,7*,F3.3,4X,FF.3,9X,F6.2,"X,F7.2,4X, 
lEe.2,4x,£8.2 ) 

GO TC 4525 
4523 SH=1. 

4525 CCNTINUL 
IF(SH)453C, 4530, 4526 

4526 PRINT 4532 

4532 FORMAT (// I X , 7HF ANGERS /) 

PRINT 4527 

45270FCRMAT I2X, 9HCCMPCNENT, 1 4 X , ^ HCLX , I 4X , 3HCL Y , 1 4X , 3HC LZ , 14X , ?HC I X , I 4 X 
1 ,3HCTY, 14X,5HCTZ /) 

4528 DC 4530 1=1 , NS 
IF(P(I 1 14530,4530,4529 

4529 PRINT 4531 I , CL X ( I ) , C LY ( I 1 , CLZ ( I ) ,CT X II ) , CT Y ( 1 ) ,C T Z ( I j 

4 531 FORMAT (3X, 1 .3, 1 5X , F 1 0 . 2 , 7 X , F 1 C. 2 , 7X ,F 1 C . 2 , 7X , F 1 0 . 2 , 7X , F 10 . 2 , 7X , F 1 C . 
12 ) 

4530 CONTINUE 
PRINT 446 

448 FORMAT (//2CH BOUNDARY CCNDITICNS 1 

PRINT 440 (SVI(I,1) 1=1,61 

449 FORMAT ( 5H SVI=6F3.0) 

PRINT 451 (SVE( 1,1) 1=1,6) 

451 FORMAT ( 5H SVE=oF3.0) 

1F(NBR)4511 ,451 1 ,4502 
4502 IF(JK)45C5, 4505, 4515 

4505 PRINT 4510( (SVIPCl, J1 1=1,6) J=2,NNN) 



OCiCCO 
0C3G1 C 
003820 
00 38 30 
0C3’J40 
0C385C 
003P60 
0C3870 
003880 
003E9C 
0C39C0 
0C391O 
003920 
00393C 
003940 
0C3950 
003960 
003970 
003980 
003990 
0C4000 
00401 C 
0C4C20 
004030 
004040 
004050 
004060 
004070 
004080 
004090 
004100 
0041 10 
004120 
004130 
004140 
004150 
004160 
0C4170 
004180 
0CU190 
004200 
004210 
0G422C 
004230 
004235 
0C4240 
0C4250 
004260 
004270 
004280 
004290 
004300 
00431C 
004320 
004330 



J 




4510 FCRVAr UF3.0 ) 004340 

GC TC 4511 OC435C 

4515 PRINT 4517( (bVOPI I* J) 1=1.6) J»2,NKN) OC4560 

4517 FrRyAT {6F3.C1 004370 

4511 CCNIINUC 0C4380 

nr 446 ^'=1,^H 0C43<50 

446 OPECPf'(M=0.0 004400 

nr 447 L=1.LM 004410 

CPtGA(U=C.O 004420 

447 nEHL)=0.0 004430 

OELC«C.O 004440 

rPECAC=0.0 004450 

IF{N'N)453.453,452 0C4460 

452 MN«0 004470 

L=1 004480 

|/=1 004490 

AA»0.0 0C45C0 

BB=O.C 004510 

CC=0.0 004520 

GG=0.0 004530 

KK»0 0C4540 

AOOT'*F.nCM 004550 

CR=BR 004560 

OMEGA! 1 )=OMtGAI 004570 

OMEGAr=OMEGAl ' 004500 

DOMsAOOM 004590 

GO TC 998 004600 

453 JJ»JJ*1 004610 

DO 454 1=1, KS 004620 

P(l)a0. 004630 

454 2( I )*C.C 004640 

GO TO 450 004650 

tND 004660 

C 004670 

C 004680 

SUBROUTINE SUBSEC { RHO , TF E T A , C , TL , HM , Z ) 004690 

RHOO«A0SF(RHOI 004700 

IF{HM-5.0)8,a ,7 004710 

7 HM«3.0 004720 

8 RA-TL/n 004730 

IF(RH0n-0. ) 3, 3,4 004740 

3 IF(RA-1. )5,5,6 004750 

5 2«0. 004760 

RETURN 0C4770 

6 lF(RA-3.) 15,15,16 004780 

15 1 * 2,0 004T90 

■ RETURN 004800 

16 lFlRA-6. ) 17, 17,18 004810 

17 2*3.0 004820 

RETURN 004830 

18 IF{HM-3. 11, 1,2 004840 

1 2*6.0 004850 

RETURN 004860 

2 2*2.0»HM 004870 
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RETURN 
4 Pn=KHCC/C 

1F(RL-1.)19,1<5,2 3 
19 1=0 
RETURN 

23 IF(RL-3. )?4,24,2S. 

24 IF( THETA-. 44)26, 26, 27 
26 2=0. 

RETURN 

27 IF(TH:TA-1.37)2F., 23,29 

28 Z=2.0 
RETURN 

29 Z=3.0 
RETURN 

25 IF(RC-6.)30,3C,i1 

30 IF(THETA-.4u)32,32,33 

32 Z=2.C 
RETURN 

33 IF(TFLTA-1.57) i4,?4,35 

34 Z=3.0 
RETURN 

35 Z=4.C 
RETURN , 

31 IF(THETA-.44)36,36,3 

36 Z=2.0 
RETURN 
END 

C 

C 

SUPRCUTINE DISTV ( AMU , A I V ,0P EGA , I L ,R , C, D I , G , SO, R I , E , FT AC , U ) 
DIMENSION U(6,6) 

AR A=( 3. 141592654/4.0) *(C*» 2-01 ••2) 

AREA=ARA«FFAC 

eE=TL-CMEGA*SCRTF(AMU/(E«ARA) ) 

P4=R«AMU»TL**4»CMEGA**2 
IFIRI ) 1, 1,2 

1 T=R»AMU*(AIY«CvrGA*TL)**2 
GO TC 3 

2 T=0.0 

3 IF(SC)4,4,5 

4 S=(AMU*(CMEGA.TL)**2)/(G»AREA) 

GO TC 6 

5 S=0.0 

6 SPT=S+T 

VV=SCRTF(P44(S-T ) •*2/4.0) 

AL1=SCRTF( VV-.5*SPT) 

AL2=SCRTF( VV4 .5*SPT ) 

CL=1 .0/1 AL1**2 + AL2**2 ) 

CHL1=CCSH(AL1 ) 

SHL1=SINH(AL1 ) 

CL2=CCSF( AL2) 

SL2=SINF( AL2) 

SPE=SINF(RE) 

CBE=CCSF( BE ) 



004330 
0C489C 
0049CC 
0C49)C 
0C4920 
004930 
004940 
004950 
004960 
0049 ro 
004 930 
004990 
0C5000 
00501c 
005020 
0C5030 
0C504C 
005050 
005060 
005070 
005C8C 
005C9C 
005100 
0051 10 
005120 
005130 
005140 
005150 
005160 
00517C 
005180 
005)90 
005200 
005210 
005220 
005230 
005240 
005250 
0C526O 
005270 
005230 
005290 
005300 
OCSilO 
00532C 
005330 
005340 
005350 
005360 
005370 
005380 
005390 
005400 
005410 
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c 

c 



CrO=CLM AL2»*:«ChU+AL U»?*CL?) 

r.Cl =CL»( ( AL^«»2»5I-L I ) / M l + ( OL 1** ‘S- 2 ) / iL. ) 

CC2=CL»(Chl 1-CL2 ) 

C03=CL<t( SHLI/flLl-SLi/ALi ) 

[.C 7 J=l,6 
l.'C 7 K=1 , A 
7 UU,K)=0.r 
L(1 , 1 I=lHE 

U( 1 ,6) = IL»SPt/(PL»A^A*l£) 

u( 2 , 2 )=coc-s*cc-: 

U(2, i)=TL*(CCl-SHT*C03l 
U(2 )=R*C02« TL»»2 

U(2,S) = (R»TL»»3/Pti|»((P‘4+S*«I)*CL3-S*CCl) 

U(3,2)=PU*CC3/TL 

L'(3, i)=CCC-T»CC2 

U(3,P )=TL*R»(C01-T»C03) 

U(3,5)^U(2,4) 

PO= 1 .C/R 

U(U,2)=Pi^*RR»CC:/TL**r 
U(4,3)=RR»( (PJ^ + r»*2 )*CC3-T*CCn/TL 
LM4,4)=U(3, 3) 

U(4,5)=L(2,3( 

U(5,2)=P4*RR» (CC 1-S*C03) /TL»*3 
U5,3)=L(4,2) 

U(5,4)=U(3,2 ) 

U(5,5)=4(2,2) 

U(6,n=-AMU«TL»CyEGA»»2»SHE/HF 

U{6,0=CRE 

RETURN 

END 



SURRCUTINE CISTVC (AKU,AIX,A]Y,OPEGA,TL,flJT,R,D,DI,G,SD, RI,EFAC,U 
DIMENSION U(6,6) 

V«=l ./(G* AJT ) 

AREA=I 3. 141 5926 = 4/4 .0 ) • ( r.«*. -L’I**2 )• FFAC 
BF = +SCR1F(AMU • ( T L » A I X *CM EG A » *» 2 • U ) 

PU=R«AMU • TL ••4*0MEGA»«- 
IFIRI >2,2,3 

2 T = R*AMU» ( A1 YtCF'EGAMTL ) »»2 
GO TC 4 

3 T=0.0 

4 IF(SC)5,5,6 

5 S*( AMU«(CMEGA«TL)»*2)/(G»ARLA) 

GO TO 7 

6 S=0.0 

7 VV=-»SCRTF(P4+(S-T)'i«2/4.C) 

SPT=s+r 

ALl= + St,RTF( VV-.5»SPT) 

AL2=+SCRTF( VV+. =bSPT ) 

CL=1.0/IAL1»»2+AL2««2) 

CHL1=CCSH(AL1 ) 

SHL1=SINH»AL 1 ) 

CL2=C0SF( AL2) 



C0542C 
005430 
005440 
00545C 
00546C 
005470 
0C54R0 
0CS49C 
0055CC 
00551C 
0C5S20 
0055 30 
00554C 
0C555C 
005d60 
005570 
0055P0 
0C5590 
005600 
005610 
005620 
005630 
005640 
0C565C 
00566C 
0C567O 
005630 
00569C 
005700 
005710 
00572C 
00573C 
0C5740 
005750 
005760 
005770 
005780 
005790 
0058CO 
005810 
005320 
005330 
005040 
005250 
OC5060 
0C5870 
005830 
005890 
0059CO 
005910 
00592C 
005930 
005940 
005950 
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SL2=SlNf ( al: ) 

SPE = SINF ( EE I 
CPE^CCSF(EE) 

CC0=CL«(AL2*»Z*CHL1 + ALl«*. «CL2) 

CC1 =CL*( ( ALL-«»2«SF L I )/AL 1 ♦ ( AL 1 • * 2« S L 2 ) / AL . ) 
CQ2=CL.(CHL1-CL2) 

C03 = CLMSHL1/AL1-SL2/AL2 I 
nr 1 j=i,d 
no 1 K= 1,6 

1 U(J,K)=C.C 
tn , 1 )=ct'E 
U(,1 ,2> = TL«W*SFE/f<E 
u(2,n=-Aru • TL •Aix**:*sPE/aL 
U(2,2»=CEt 
U(3,i)=CCC-S*C02 
U(3,4)= TL •(CCl-SF'r»COi ) 

U( ?,5)=R*CC2«TL«»2 

U(3,6) = (R*TL»*?/P'>)*(-S»CC1+(P4 + ^«*2)*CC3) 

U(U,3)=PU*CCa/'TL 

U(U,4>=CnC-T«CCL 

U(U,5)= TL»R*(CC1-T*C03) 

LM4,6)=0(3,5) 

RR=1 .C/R 

U( 5* ?)=P4 »RR*CC2/TL»»2 
U(5,4)=RR»(-T*CCl+(P4+r**2)«CO?)/TL 
U ( 5 , S ) =L ( 4 , 4 ) 

W(5,6»=U(3,'*) 

U(6,,3)=P4*RR» (cr 1-S*cn3 )/TL**3 

U(6,'4 )=U(S, 3) 

i;(6,i)=L(4,i) 

U(6,6)=U(3, i t 

RETURN 

END 



SUBRCLTINE SEIELC ( f I , TL , P , / , G , SI , 0, E , FT Af. , A ) 
OINlNSIGN A(6,6) 

AR A=( 3.14 1502654 /4.C) M E*« 2-E. 1 ) 

AREA=ARA.FFAC 
UL =TL/6 
DO 1 J=2,5 

1 A(1,J)=0.0 
A ( 1 , 1 ) = 1 . 0 

All ,6)=0 L/( l*ARA) 

A( 2 , n=o.o 

A«2,2)=1.C 

A(2,.‘)=UL 

A(2,4)=DL*«./«R/2.r 

IF(SU)2,2,3 

2 A(2,5)=LL«*A»R/6.C-DL/'«G*ARLA) 

GO 1C 4 

3 A’(2,5)=UL*»r,«R/'6.r 

4 A(2,61=C.O 
A(3, 1 )=C.O 



0C5960 
005970 
0059EC 
005990 
0C60C0 
006010 
0C6Ci.C 
006030 
0C60M0 
0C6Cd 0 
00606C 
006070 
006080 
00609C 
006100 
0061 10 
006120 
006130 
006UC 
006150 
0C6160 
006170 
0061SC 
006190 
006200 
006210 
00622C 
006230 
006240 
006250 
006260 
0C6,i70 
006280 
00629C 
006300 
006310 
006320 
006330 
006340 
006350 
006360 
00637C 
006380 
006390 
006400 
006410 
006420 
006430 
006440 
006450 
006460 
006470 
006480 
006490 
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I 







A( i,2 )=C.O 
A( ?, 3 ) = K0 
A(?,4)=uL»R 
A(3,5)=A(r,L) 

A(i,6 )=0.0 
Dr; 5 j=l,e 
A(U,J)=C.C 
A(5, J)=O.C 
5 A(6, J)=O.C 
A ( U , 4 ) = 1 . C 
A(4,5)=(iL 
A(5,5) = l .C 
A(6,6)=1.0 
RLTDRN 
bND 
C 
C 

SURKCUriNE SFIELC ( D I , I L , A J I , R , Z , G ,S UR 1 , 0, FF AC , A ) 
DIMLNSICN A(6,6) 

AR EA=( 3. 14 1592654/4.0) •( r.«*;-r I »»2 )• FF AC 
UL=TL/Z 
DC 2 J=1 ,6 
r,G ^ K=),6 
2 A(J,K)=0.0 
AM, 1) = 1.0 
A ( 2 , 2 ) = 1 . C 
A(3,3)=1.0 
A14,4 )=1 .C 
A(5,5)=1.0 
A(6,6) = l .0 
IF(5CRI-0. )7,7,R 

7 A(3,6)=DL«*3«R/6.C-CL/(G»ARtA) 

GC TC 1 

8 A(3,6)=nL««3*R/6.C 
1 A( 1 ,2)=UL/IC*AJT) 

A(3,4)=r.L 

A(3,5)=DL**2»R/2.C 

A(4,5)=UL»R 

A(4,6)=A(3,5) 

A(5,6)=(.L 

RETURU 

END 

C 

c 

SUPRCUTINl CFIELU (R,Z,G,SD,RHC,THErA,n,OI,E,A) 
DIMENSION A(6,6) 

PHI=THETA/Z 

RD=AeSF(RHOI 

1F(RHD1,1,2 

1 V=-1.0 
GO TC 3 

2 V=1.C 

3 CP=CCSF(PFI ) 

SP=SINF(PHI ) 



C06500 
00651 C 
0 C6 j C 

cc6d;c 

0C654C 
00655C 
0C6560 
OC6570 
006580 
00659C 
006600 
006610 
006620 
0C663C 
006640 
0066 jO 
006660 
006670 
OC66‘»0 
OC6690 
006700 
006710 
006720 
006730 
006740 
006750 
006760 
006770 
006780 
006790 
OC68CO 
006810 
006820 
006830 
006840 
0C6850 
006C6C 
006870 
006830 
006890 
006900 
006910 
0C6920 
006930 
006940 
006950 
006960 
006970 
00698C 
OC6990 
007000 
007010 
007C20 
OC7C30 
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1 r 4 j=i ,6 

cr. u K.= i,6 

4 A(J,K)=C . C 
A( 1 , 1 )=CP 
A( 1 ,2)=-V*SH 
A( 1 1 .0-cr*) 

A(1 ,4) = -v«Ri;««2«r<»(PHi-SM 

AP= <. 14 1=924 b4«(L«».?-C I )/4.C 

U=1 .C/(c«A'rO 

F?= I SP-I'HI»cP )• .b 

Fl = (PhleCP + ‘^P)*.5 

F5=l2.C-2.0»CP-i^l-I*SP)*.‘; 

r6=(2.0«PHH-PHI«CP-i.0*SPl».5 

An, 5 )=v*(Rr;*u»PHi«sp».:-R*F 5 » 9 r»» 

A(1 ,6)=RD«*V*Fl + ;l«F6*'<n»<i 

A(2, 1 1=V»SP 

A(2,£)=CP 

A(2,3)=KD«SP 

A(2,4)=R«R0»»2»(1.0-CP) 

A(2,5) = :U;*F'»(vn + r«rc*»2) 
Al2,6)=An, 

A(?,i)=1.0 
AI3,4 )=RO«R«PFI 
A(3,5)=A{2,4| 

A(?,6)=A( 1,4) 

A(4,4 ) = 1 .0 
A(4,5)=M2, 7) 

A(4,6)=An, A) 

A(5,5)=LP 
A(5,6)=A( 1,7) 

A(6,5)=A(2, 1 ) 

A(6,6)=CP 

RGTLRN 

END 



SUBKCUTINL CFIELC ( A J T , P , Z , C , SCR I , 
DlMtNSICN A(6,5) 

PHI = T)iCT A/Z 
RhCD^RHC' 

IF( RHC-C.O) 11,11,12 

11 V=-1.C 
GO TC 1 

12 V=1 .C 

1 CP=CCSF(PHI ) 

SP=SINF( PHI ) 

U=1 ,0/( C *AJT ) 

Fl= ( PHI«CP + SP )/2 .C 
F3=(SP-CP*PHI )/?.(' 
F5=(2.C-2.0»CP-PHI«SP)/2.C 
F6=(2.0«PHl+PFI»CP-3.r»SP)/2.C 
RHO=ARSr (RHC) 

, DO 2 J=1 ,6 
DO 2 K=1 ,6 



007C40 
0C7C5C 
COTCtO 
007C70 
n07C- C 
0C7C9C 
007 ICO 
0C7’. I 0 
007 12C 
007170 
C07140 
0C7):;C 
0C7if^C 
007170 
0071 CO 
00719C 
0C7cC0 
007' 10 
007220 
007^50 
007240 
0072S0 
007240 
007^70 
007230 
CC7290 
0O7S00 
007.'51C 
OC73"0 
007ii0 
007740 
007iSC 
007360 
007i70 
007380 
DC739C 

, THETA,D,DI,FF AC, A) 007400 

007410 
007420 
0074i0 
007440 
007430 
0C746C 
007470 
007430 
007490 
007300 
007S10 
0C7520 
00733C 
007540 
007550 
0C7560 
007570 
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2 A(J,K)=C.C 
If (bCRUl’, 

13 ARE,A={ 3. 141b9 265U/U.C)*(C«*:-DI*i!2 )» ff AC 
Sn=RHC« THETA/ ( G»Z»AREA) 

GO TC lb 

14 SL)=O.C 

15 A( 1 , 1 )=CP 

A ( 1 ,2) = ( K*RHC;»F1 )-(.RHO»F3«R ) 

AM ,4)=V«SP 

A(1,5)=V*(W+R)»(RHO*RHI»SP1/1.0 
A( 1 ,6) = V« (W + R ) •( RHC**2)»F? 

A(2,2l=CP 
A(2,5) = \/«SP 
A(2,6)=V*RHC*( 1 ,0-CP) 

A(3,11--A(2,6) 

A(3,2I=-1.0«A( 1,6) 

A(3,6)=1.C 

A(3,U)=RHC*bP 

A(3,5)=((R*(RH0»»2)«PH1*SP1/2.C)-W*RHC**2»F5 
A( 3,6) = ( R*RHO**3«F3)-U«F6*RHC** ^)-SO 
A(4 , 1 )=-y/»SP 

A(4,2)=-V*(t*+B)«RH0»PHI«SP/.'.0 

A(4,4)=CP 

A(4,5) = (R*RHC»F1 )-( In»RHC«F3 ) 

A(4 ,6>=A(3,5) 

A15,2)=-V*SP 

A(5,5)=CP 

A15,6)=A(3,4) 

A(6,6) = l .C 
RHO-RHCC 
RETURN 
END 
C 
C 

SUBROUTINE POINT ( A^U , A I V , OR EGA , TL ,Z , R I , B ) 
DIMENSION B(6,6) 

AM=AMU*TL/( Z-1. ) 

UC 1 J=1,6 
DC 1 K=l,6 

1 BU,K)=K).0 
B( 1 , 1 ) = 1 .0 
P(2,2)=1.C 
B(3,5) = 1 .C 
IF1RI )2,2,3 

2 R(4,3)=-AM*(AIY.CyEGA).*2 

3 8(4,4)=1 .0 
B(5,2)=AM*0MFGA*«2 
B(5,5) = 1 .0 

P(6, 1 )=-B(5,2) 

B(6,6 ) = 1 .0 
RETURN 
END 
C 

c 



0C75RC 
0CZE90 
0076CC 
0076 1C 
0C762C 
C07tiC 
00764C 
0C76‘;0 
0C7660 
00767C 
00763C 
00769C 
0077CC 
0C771C 
00772C 
0C775C 
Q0774C 
0C7750 
007760 
0C7770 
007780 
007790 
007800 
0C7810 
007820 
0C7830 
007840 
0C7850 
007060 
007870 
007880 
007890 
007900 
007910 
007920 
007930 
0C7940 
007930 
007960 
007970 
007980 
007990 
008CC0 
008010 
008020 
008030 
008040 
008050 
008060 
008070 
008080 
008090 
008100 
008110 
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SUBRTUTINf PCINC ( A VU , a I X , A I Y , OP ■ G A, T L , ^ , SD I , P I 
niMENSILN P (fc,6) 

AP=APU»TL/(Z-1.C) 

DO 1 J=l,6 
DO 1 K=1,t 
1 P(J,K)=C.C 
P ( 1 , 1 ) = 1 .C 

B(2, n=-AP*( AIX»n«ECA(»»Z 

B(2,2)=1.0 

P ( 3 , i ) = 1 . D 

B<4,A)=1.0 

lF(SDRI-0. )9,9, IC' 

9 Pt5,4)=-AP»(AIY«CNEGA)«»2 
GO TC 11 
10 b(5,A(=0 
n EI5,£)=1.C 

P(6,3)=AM»CP[:CAii* 2 
P(6,6)=l .0 
RETLRN 
END 
C 
C 

SUBROUTINE PATMUL (A,B,Z,U) 

DIMENSION A(6,61,P(6,61,C(6,6),D(6,S),U(6,61 
DC 1 J=1,6 
DO ) K=1,Z 
C ( J,KI=O.C 
DO 1 L=l,6 

1 C(J,K l=U( J,K)+8IJ,L )*A(L,K) 

N=Z-1 .8 

DO 2 J=l,6 
DO 2 K=l,6 

2 U(J,KI=C(J,K) 

IF(Z-2.0)6,6, 3 

3 DO 5 M=1,N 
DO 4 J=1 ,6 
DO 4 K=l,6 
D(J,K)=0,0 
DO 4 L*1,6 

4 0(J,K)=D(J,K)+CIJ,L)*U(L|K) 

DO 5 J=1 ,6 

DC 3 K=l,6 

5 U(J,K)=D(J,K) 

6 DO 7 J=1 ,6 
DO 7 K=1 ,6 
D( J,K)=O.C 
DO 7 L=l,6 

7 C(J,K1=D(J,K)+A( J,L1*U(L,K) 

DC 8 J=1 ,6 

DO D K=1 ,6 

8 U(J,K)=D(J,K) 

RETURN 

END 



0C8 1 20 

ccei 3c 

0C81U0 
0C31 bO 
0C8U0 
0081 fC 
0C8180 
0C819O 
0082C0 
OC82 10 
0C8220 
008230 
008^40 
0C82S0 
0C82O0 
008270 
008^80 
008290 
008300 
008,310 
008320 
008330 
OC034O 
008350 
008360 
008370 
008380 
008390 
008400 
008410 
008420 
008430 
008440 
008450 
008460 
0C8470 
0084SC 
008490 
008500 
008510 
008520 
008530 
008540 
00855C 
008560 
008570 
008580 
008590 
008600 
008610 
008620 
00863C 
008640 
008650 



D-17 



SL.'BKCurr>iE ?icn; TL,cyuGA,A iv,-u ,l ) 

DIMENSION U(6*6) 

AM= AMJ< TL 
no 1 J=],6 
DC 1 K=l,6 

1 D{J,K)=D'.C 
U( 1 , 1 ) = 1 .0 

0(6,1 l=-flN»r»'LG4»»? 

U(6,6)=1.C 

0 ( 2 , c ) = 1 . 0 

0(2, J)-IL 

0 ( 3 , i ) = 1 . C 

0(H,2)=AM«TL»CM£GA»«2/2.r 

1 F ( R I ) 2 , 2 , 3 

2 0(U,i)=AM.Ci^eCA»»2i.(TL»»'/t.0-AIY».2 ) 

GC Tt h 

3 0(4,3)=AM.(CMEGA.TL)»*2/6.C 

4 0(4,4 ) = 1 .0 
0(4,5)=IL 
0(5,2)=-0(6, 1 ) 

0(5,31=0(4,2) 

0(5,51=1.0 

RETURN 

LND 

C 

C 

SOPRCOTINb RIGIC ( AMU, TO , A I X , OMEGA ,A I Y , A JT, SCR I , 01 
DIMLNSION 0(6,6) 

DC u. J=1,6 
DO 2 K=1 ,6 
2 U(J,K 1=0.0 

am=amo*tl 

0(1,11=1.0 

0(2,1)=- AMU •TL*(AIX«CME3A)*»2 
0(2,21=1.0 
0(3,31=1.0 
0(3,4 )=T0 
0(4,41=1.0 

0(5,3 l = AMoTL »CM EGA** 2/2 .C 
IF(SDRI-C.) 30,30,31 

30 U(5,4)=AM*OMEGA**2*(TO**2/6.0-AIY**2 1 
GC 1C 1 

31 U(5,41=AM*CMEGA**2*(IL**£/6.0) 

1 0(5,51=1.0 

0(5,6 )=T0 

0(6,3 1=AM*0MEGA«*2 
0(6,41=0(5,3) 

0(6,61=1 .0 

RETURN 

END 



SOBRCUTINE STIFCC ( THET A , RFC , 0 1 



.0C6C60 
0C8670 
008630 
0086 ')C 
0C87CC 
003710 
0087:0 
0Q875C 
OC87-.0 
OCS75G 
008760 
0087 70 
008730 
008790 
OCS5CO 
008- 1 0 
0088^:0 
0088^0 
0C3J40 
00BH5C 
0CB56C 
008670 
008680 
008890 
0C89C0 
008910 
0C892C 
008930 
008940 
008950 
006960 
OC6970 
008930 
008990 
00900C 
009010 
0C902C 
0C9C3C 
009040 
009050 
009060 
009070 
009C8C 
009C9C 
009100 
0C911C 
009120 
009130 
009140 
0C915C 
009160 
009170 
009180 
009190 
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i 



D1ML\S1C\ U(6f6) 
f.r I J=1,6 
rr ! K=i,6 

1 l'{J,K)=O.C 

2 V=-1.C 
GC 1C, 

5 v =].0 

4 CT = CCSF( THlCA ) 

ST=SINF( TFETA ) 

DM , n=GT 
Li( 1 ,2 )=-V»ST 
U(2, n=V»ST 
L(2,2 )=CT 
U(3,i)=l .C 
IJ ( 4 , 4 ) = 1 . C 
L'(5,5)=LT 
L(5,6)=-V*ST 
U(6,E )=V*ST 
U(6,d )=CT 
RETURN 
END 
C 
C 

SUPRCUTiNE STIFCO ( TH ET A ,RhC , U ) 
CIMLNSICN U(fc,6) 

IF(RhC) 1,1,2 

1 V = -I.O 
GO rc =1 

2 V=1.C 

3 CT=CCSr(THLTA ) 

ST=bINF(THETA) 

DO 4 J=1 ,6 

CO 4 K=l,6 

4 U(J,K)=O.C 
UM , n=CT 
U( 1 ,4 )=5T*V 
U(2,2)=CT 
U(2,5)=ST«V 
U(3,?)=1.C 
L'(4,1)=-ST*V 
U(4,41=Cr 
U(5,2)=-ST»V 
U(5,5)=CT 
U(6,6) = l .0 
RETURN 

END 

C 

C 

c 

SUBRCLTINE HANGER { CL X , C L Y , GT Z , U 
DIMENSION U{fc,6) 

DO 1 J=1,6 
DC 1 K=l,6 



org^ico 
OCR; 10 
0C92:. c 
OCR.: ^0 
OCR^UO 
OCR. -■€ 
OC^OlcC 
oor: cc 

0CR2 -0 
0CR2 RC 
OCRICC 
OCR^IC 
OOR ^20 
OCR:.?C 
0CRa4C 
OCR^'iO 
0CR;60 
CCRf ’0 
OCRiSC 
OCR.-'RO 
0CR400 
0CR41C 
0CR42C 
0CR4 iO 
0CR44C 
0C945C 
OOR 4 oO 
0CR470 
0CR4RC 
OORURO 
OCRbCO 
0CR510 
OCR-.ZO 
009230 
OCR540 
0095SO 
OCRbdC 
00957C 
OCRbRO 
OORbRO 
00R6CC 
OORolO 
0C962C 
009650 
0C964C 
009650 
0CR66C 
009670 
OOR6CO 
0C96RC 
0C9700 
009710 
OOR72C 
00R73C 
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c 

c 



c 

c 



1 L(J,K)=0.0 

un , 1 ) = i .r 
u ( 2 , 2 ) = ’ . 0 
U ( 3 , £ ) = 1 . C 
U(«j , 3)=CT2 
U(i» ,4) = 1 .C 
U(5,2)=-CLY 
U(5,5)=1.C 
U(6, I 1=CLX 
U(6,6)=l .0 
RETUi^N 
END 



SUB^CUTINE HANGlC ( C T X , C L Z , C TY , L ) 
DIMLNSICN U(6,6) 

DC 1 J=l,t 
DC 1 K=1 ,6 
1 b(J,K)=0.C 
D < 1 , 1 ) = 1 . C 
U(2, 1 )=CTX 
U(2,2) = l .0 
U ( 3 , S I = 1 . 0 
i; ( u , 4 ) = 1 . c 
U(5,4I=CTY 
U(5, 51 = 1.0 
U(6,5 )=-CLZ 
L(6,6)=l .0 
RETURN 
END 



SUPRCUTINE PRANCH ( SV P , VV , PE I , I ,U 1 
0DIMcNSICNSVB(6,22), VV( 6,6),U(6,6),R1( 
1 R ( 3 , 3 ) , G 1 ( 3 , 3 ) , C 2 ( ? , £ ) 

P=1 

DC i N=1 ,6 
IF(SVetN,M 1 ) £,£, 1 

1 DO ^ L=1,6 

2 VV( L,r') = VV(L,N) 

M = y+1 

3 CCNTINUc 
DO 4 L = 1 ,3 

R1 ( 1 ,L)=VV( 1 , L) 

P1(2,L)=VV(2,L) 

R1 <5,L»=VV( 3,L) 

R2( 1 ,L)=VV(4,L) 

R2(2,L)=VV(5, L) 

4 R2(3,L)=VV(6,L) 

CALL INVERT ( R 1 , 3 , D , LL , 

DC ^ J=l,3 
DC K=l,3 
R ( J,K)=C.C 
DO 5 L=1,3 



009740 
C097SO 
009760 
0C9770 
0C9 r?c 
0C979C 
OC90CO 
0C9D10 
009t^C 
009830 
0C9L'40 
0C9f. £C 
009-60 
009E7C 



009E9C 
00990C 
00991C 
009910 
0099 30 
00994C 
0C9950 
0C9960 
009970 
0099=^0 
009990 
010000 
010CIC 
01CC20 
010C30 
0 10C40 
01C05C 
010C60 

25) ,LL (25),NP(25 I ,R2(3,3 ), 010070 

010080 
01 0090 
010100 
010110 
01012C 
010150 
010140 
010150 
010160 
0101 70 
01018C 
010190 
010200 
010210 
010220 
010230 
010240 
01025C 
010260 
010270 
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1 



5 R(J,K)=r(J,K)+R2(J,L)*^l(L|K) 
SP=i,INF( PMI ) 


.010250 


CP=CCSF ( PHI J 


010290 


Cl n, 1 )=CP 


010300 


Gi ( 1 )=-sp 


01 0.1 1C 


Cl ( 1 ,i ) = C. 


010320 


Gi(i,n = sp 


0 lO.'.iO 


G1 (^,l)=CP 


01C34C 


Gl(i,1)=0. 


0103-C 


Gi (i, n=o. 


C1036C 


Cl (3,21=0. 


010370 


G1(3,3)=1. 


0103E0 


G2( 1 , 1 ) = 1 . 


010390 


G2( 1,2)=0. 


OlOi4CO 


G2( 1 ,3)=0. 


010410 


G 2 ( 2 ,n= 0 . 


01 04"0 


G2(2,2)=CP 


010430 


G2(2, <) = -SP 


010440 


G2( /. ,1)=0. 


010450 


G2( 5,2) =SP 


010460 


G2(i,i)=CP 


010470 


DO (, J=1 , 3 


010480 


DO t K=l,,3 


010490 


R 1 ( J , K ) = C . 0 


010500 


DO 0 L= 1 , 3 


01 05 10 


> RHJ»K1=R1(j,k)+C2( J,L1«P(L,KJ 


010520 


DC 7 J=l,3 


010530 


DP 7 K=1 , 3 


010540 


R( J,K)=C.O 


010550 


DO 7 L=l,3 


010560 


R(J,K»=R(J,k)+R1(j,l)*GI (L,K) 


010570 


DO £ J=l,t 


0105PC 


DO e K=l,t 


010590 


0( J,K)=C.C 


010600 


L) n , 1 ) = 1 . c 


010610 


0(2,21=1.0 


010620 


0(3,31=1.0 


010630 


0(4,41=1.0 


010640 


0(5,51=1.0 


010650 


0(6,61=1.0 


010660 


DO 9 L=1 , 3 


010670 


0(4,L1=R( 1,L1 


010680 


0(5,L1=R(2,L1 


010690 


U(6,L1=R( 3,L 1 


0107C0 


RETURN 


010710 


END 


010720 


SUBROUTINE STflVEC ( SV, N, 8C 1 


010730 

010740 

010750 


DIMlNSICN BC(6,221 


010760 


DO 10 J=l,6 


010770 


BC( J,N1=0. 


010780 


IF( SV-1 . 1 1 , 1 ,2 


010790 


BC(4,M=1.0 


010800 




oioeic 
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c 

c 



per, ,\) = i.o 

eC(6,M = 1.0 

kcruPN 

2 IF(iV-2. )■', ;,4 

3 PC ( 1 , M = 1 . 0 
RC(.- ,N) = 1 .0 
FC( ;i,N)= 1 .0 
PLTLKN 

4 IF(SV-3. )5,t>,6 

5 PC(.->,N) = 1 .C 
PCr.,N) = l .0 
EC(C>,M = 1.0 
RETURN 

6 IF(bV-4. 17, 7, F 

7 P,C( 1 ,M = 1 .0 
hC( : ,.\) = 1 .0 
PC ( b , N ) = 1 . 0 
RETURN 

8 PCC: ,N)-1 .0 
PX ( : , N ) = 1 . 0 
PC(6,N)=1.0 
RETURN 

END 



SIJERE'UTINE BRANCC ( SV f' , VV , P E I , I , U I 
GDI PENS iriNSVH(6,22) ,U(6,t»,VV(6,6»,V( t,i),Rl(25,:5l,R2(3, 3),R(3o), 
1G1 (3,0) ,U2( 5, 3) ,S(3,3 1 ,LL12b),V’^(25) , SR(6,2?),UV(6,o) 

P=1 

SR( 1,PI)=SVF ( 1,RI) 

SR(2,P1)=SV'M3,M) 

SR( 5,M )^SVE(4,’'I ) 

SR(4,P1 )=SVE(2,y I ) 

SR( 3,PI )=SVL-( 5,M ) 

SR(&,M)=SVP(6,PI) 

UV( 1, 1 )=VV( 1,1) 

UV( 1 ,2 )-=VV( 1,3) 

UVt 1 ,3)=VV1 1,4) 

UV(1,4)=VV( 1,2) 

UV( ) ,5)=VV1 1 ,5) 

UVl 1 ,6)=VV(1 ,6) 

UV(^, 1 ) = VV( 3, 1 ) 

UV ( , : ) = VV ( 3 , 2 ) 

UV(2,3)=VV( 3,4) 

UV12,4)=VV( 2,2) 

UV(^,‘i) = VV(3,5) 

UV(E,6)=VV(3,£) 

UV( 3 , 1 ) =VV( 4 , 1 1 
UVl 3,2 )=VV( 4, 3 ) 

UVt 3,^)=VV(4,4) 

UV(3,4)=VV(4,2) 

UV(o,5)=VV(4,5) 

UVt 3,6)=VV(4,£) 

UVt 4, 1 )=VV(2, 1 ) 



oio^r'o 
U1CE3C 
0 1 Q J 4 C 
u'1oe:o 

OlOl'bO 

nio!.?c 

0 ICi'PO 
010390 
010900 

oioo:o 

010920 

010930 

010940 

01095C 

0109:20 

010970 

010900 

01099C 

0 1 ICCO 
OllClO 

01 1C2C 
011030 

0 11040 

01 1090 
Cl ICtO 

0 1 1070 
011030 
011090 
011100 

0 1 n 10 
011120 
011130 
01 Hue 
011150 
011160 
011170 

01 n EC 

011190 
011200 
01 1210 

0 IIE.C 
011230 
011240 
011250 
011260 
011270 
011230 
01 1290 
011300 
01131C 
011320 

01 1330 
011340 
011350 
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f 







uv(m , j I =vv I : , I 

UV( U , ^ )r vvi ?, u ) 

uv{u,u)=vv(:,:i 

i:V(*4,: (-VV(t ,i) 

UV{ , 1 1 =vv( _ , t ) 

UV( D , 1 1 - VV( D, 1 1 
UV( J,2 )=VV( 5, i) 

UV(-;, ■>) = VV(5,H) 

uv(o,u)^vv 

UV(S, t)=VV(l,5) 

UV(ry,6) = VV(h,6) 

UV(f , 1 l=VV(fc, 1 ) 

UV(C,2) = VV(<‘,2) 

0V(C, =') = VV(6,U) 

UV(|'.,14 )=VV(6,2) 

UV(',,S)=VV«(, 5) 
UV(6,6)=VV(e,6) 

DC i N=l,6 
IF(SR(N,f’I) ) i,3, 1 

1 DO 2 L=1 

2 UV(L,^')=lJV(L,^l 
M=M+ 1 

3 CONTINUl 
no u K=r,3 

R1 ( 1 ,K)=lJVn,K) , 

R1 (i,K)=UV(2,K) 

R1 (5tK)=UV( 3,K) 

R2( l,K)=UV(u,K) 

R2(2,KJ=UV( 5,K) 

4 R2(:,,K)=UV(6,K) 

CALL INVERT ( S 1 , 3 , D , L L , MN ) 
nr j j=i ,3 
CO L K=1 ,3 
R( J,K)=C.O 
DC 5 L=1t5 

5 R(J,K1=R(J,K)4R2(J,L)*R1 (L,K) 
SP=SINF( PHI ) 

CP=LCSF( PHI ) 

G1 ( 1 , 1 )=f.P 
G1 ( 1 ,21=0. 

G1 ( 1 ,3)=-SP 

Gl( 2 ,n= 0 . 

G1(2,2)=l. 

G1 (2,3)=P. 

G1 (3,1 ) = SP 
G1 (3,Z1 = 0. 

G1 (3,3)=CP 
G2I 1,1)=CP 
C2I 1,2) = SP 
G2( 1 ,3)=0. 

G2(2, 1 J=-SP 

G2(^,2)=CP 

G2(2,3)=0. 

G2(3, 1 )=0. 



'>11560 
)1 1 i70 

1)1 1 i ’0 

01 1 ;0C 
OIKCO 

0 11 >( ’ c 

01 I'^.O 
011450 

0 1 1 4 q C 

01 luEiC 
0 11460 
01 1470 
01 lu"C 
011400 
oilbOO 
011610 
011520 
011530 
01 1i40 
01 1 55C 
011=60 
011=70 
01 15R0 
011500 
01 16C0 
011610 
01 162C 
011630 
011640 
011650 
011660 
01 lc7C 
OllofiC 

0 11690 
011700 
011710 
011720 

01 17.30 
01 1 740 

0 1 1 7=0 

01 17o0 
011770 
011780 
01 179C 
011800 
01 IdlO 
011820 
011830 
011840 
011830 
011360 
011870 
01 18PQ 
011690 
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G2( ;,2) = 0. 

G2( 5,; )=1. 
r G (. J=1 ,I 
ri: o K = i ,3 
R 1 ( J , K ) = C . 0 
DO f. L=l,? 

6 Rt(J,K)^Rl(J,K)+C:(J,L)»R(L,K) 
Uf; /' J=1,3 

H’r 7 K=1 , 3 
5i(J,K)=L.C 
LC r L=1 , 3 

7 S ( J ,K)=S( J, <) iR U J,U *C 1 (L,< ) 

or j = i,t 
no c K = ) , ^ 

3 u ( J , K ) = 0 . 0 
U( 1 , 1 ) = 1 .0 
U(2. 1 )=3( 1, n 
L ( 2 , 2 ) = 1 . 0 
0(2,51=311,2) 

U12,U)=S(1, 3) 

U ( 3 , * ) = 1 . C 
U(U ,i4 1 = 1 .C 

ij(G, n=sc , n 

U15,31 = 5C,2) 

L'(5,4) = S(2, i) 

U(5,bl=1.0 
U(6, n=b( 3, 1 ) 

U(6,5)=S(3,?I 
U(6,li)=S(3,3l 
L ( 6 , 6 ) = 1 . 0 
RtTURN 
EMO 
C 
C 

SUPRCUTINc ST^VCC (SV,N,PC1 
OIMfcNSlON BC(6,221 
no 10 L=1,6 
10 6C(L,N1=C.0 

IF(SV-1.0)1,1,2 

1 EC(2,N)=1,0 
PC(5,N)=1 .0 
RC(6,M = 1.0 
RETURN 

2 IF(SV-2,0)3,3.i» 

3 BC( 1 ,N) = 1 .0 

eC(j,M = 1.0 

hC(U,N)=1 .0 
RETURN 

U IF(i.V-3.C15,5,t 
5 BC(^,.N 1=1.0 
eC(<4,N) = 1 .0 
BC(6,N)=1 .0 
RETURN 



0 1 1')G0 
011910 

01 1 R'C 
01 1930 
0 1 19U0 
0 1 19-jO 
01 1 9iC 
0 1 1970 
011930 
C 1 1990 
012CCC 
012C10 
012020 
012C30 
012C4C 
012CSC 
0 12C60 
012C70 

0 12C :c 
012C9C 
012100 
012110 
012120 
012150 
0 1214C 
012130 
0 1216C 
012170 
012180 
012190 
0 122CC 
012210 
012220 
012230 
0122U0 
012250 
012260 
0122 70 
012280 
012290 
012300 
01251C 
012320 
012550 
012340 
012350 
012360 
012370 
01233C 
012390 
0 1240C 
012410 
012420 
012430 
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7 ec( I , o) = i.i 

PCI •< , M - 1 .0 

r c ( o , 1 . c 
RC rt'- . 

C'u: 

c 

c 

SUBiiCl.'TINL li'VERT (A.m.C.L.N) 

c pRnJ;<n^■ fcr FI^CI^G Th : inv rsc f ^ Nxf.‘ matrix 
c 

D I Nc NM C r, A ( 2 E , ; t> ) , L ( 2 F ) , ■“ ( - 3 I 

C SEAiCH FTR LARGEST tLLVEAT 

tJ=1 .G 

Cl'Pt. K = l,\ 

L1K)=K 
''(K) =K 
e]G«=a(K,h 1 
(;r?0 l-K ,G. 

r.r2c J=K,N 

IF(AHSF(EIGA)-ap'^F(A(I,jni 10, 2., 20 
10 PIGA=A(l,J) 

L(K)=l 
M(K)=J • 

20 CONTINLIL 

C INT' HC^la^r,E RCk^S ' 

J = L (K I 

IF(L(K)-F) ::-,25,-2 
25 C030 1 = 1, r, 

HCLL = -a(K,1 ) 

A ( K , I I = -. ( J , 1 1 
30 A( J, I )=HCLU 
C INTlRCHANGE CCLLyEf. 

35 I=M1K1 

inf'lKl-KI 1*5, UE, ’7 
37 0040 j = l,rN 
H0LF. = -A( J,K ) 

A(J,K)=A(J, 1) 

40 A(J,I)=hnL0 

C DIVIOL COLliyN ev MNI,S PIVCT 

45 L'05 3 I = 1,N 

46 IFn-X)5C,55,50 

50 A( I ,K )=A( I,K ) /(- A(K,K I I 

55 CONTINU.; 

C REOGCL FATRIX 

006:- I = l,f 
LC65 J=1,N 

56 IF( l-K) 57,65,57 

57 IF(J-K) 60,65,60 

60 A( I , J)=A( I,K) *A(K, J 1 + M I , J ) 

65 CONTP-;L'L 

C DIVICE ROi„ hY pivot 

[.'C7h J=1,N 
68 IF ( J-K) .’0,75, 7C 
70 AIK, J)=A(K,J)/A(K,K1 



012440 
0 12** 50 
212.*6C 
012A70 
OUA'C 
0 124 0C 
0 12.-00 
012-10 



0 1 2 :• 5 c 

01254C 

012550 

0 12-J-60 

0 123 i-C 
0 125.'’C 
012500 
0126C0 
C 12o 10 
0 126 2C 
012630 
012640 
012650 
012660 
0 12o70 
(0 1 26 50 
012600 
0127CC 
012710 
0127'C 
012730 
012740 
0127^0 
012760 
012770 
012780 
012700 

oi2rcc 

012B10 
0 1 2320 
012-30 
012:40 
012050 
012o60 
012F70 
0 12e'’C 
012800 
012900 
0 1291 C 
01292C 
0 12930 
012940 
012950 
012960 
012970 
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75 Cf M ira.L 

c LrrniMj.:, pr.( tur.T Pivrrs 

nMMA(K,Ki 

L P/ PLACf- PI VCT PV '-'wC. I 

A ( K , K ) ^ 1 . t , K ) 

'iO CTN1 IMJ;. 

C riNAL TCW A'JL CCLt p" I\7EPtPANCL 

K = N 

1 00 K= ( K-1 1 

ir(Ki i:o, 150, ic.^ 

1C? I=L(K1 

ir(i-K) 1?0,120,1(5 
los urnc j=i,.n 

HOLl =Al J, K ) 
ft ( J , K ) =- A U , I 1 
lie A( J, I l=lirLU 
120 J=K(K) 

rF(J-K) 100,100,125 
125 001 3C 1 = 1 ,fti 
HflU=-'.(K , I ) 
ft(K, n=-ft{ J, 1 ) 
liO AU, I l=HOLn 
CO TC ICO 
150 Pf-.Tljn\ 

LM. 

c 

c 

su^-KouT r'f ^l^yftr nj,i.:u,i,vi 
OIMlNSIOO U(t,61,LU(c,6),V(£,61 
IF( I-1)c1 ,6 ,61 
62 l/C 70 J= 1 ,6 
re 7 0 K.= l,6 

70 V< J,K)=b( J,K) 

RETORM 

61 on 71 J=1,f; 
t;r 7 1 K=i,6 
V( J,K1=C.C 
DO 71 L=1,6 

71 V{ J,K1=\, ( J,K)+U( J,L)«UU(L,K) 
RETURN 

EMU 

C 

C 

SUBROUTINE FIMBRA { U , VV , ft ML’ R , V ) 
DIMENSION UI6.61 ,VVI6,6) ,V(6,6) 
IFI NNHR-1 11,1,3 

1 DC : J=1 ,6 
EC 2 R=1 ,t 

2 V(J,K)=U(J,M 
RETURN 

3 EC 4 J=1 ,6 
00 4 K= 1 , 6 
VIJ,K)=C .C 
DO 4 L = 1 ,6 



0 125^ C 
C 12900 

';i3cco 

C 13010 

c 1 3c:c 
Cl 30 ^0 
Cl 30 4 0 

01 3C5C 
0 1 3C60 
C13u70 
0 130-!C 
0 1 3000 
01 31C0 
013110 
0 1 31 20 
C 1 31 jO 
013140 
01 3150 
013160 
0131 ■’C 
01 31 30 
013100 
013200 
013210 
013220 
013230 
013240 

0 13.. 50 
013^60 
013^70 
cn.LEc 
013290 
013300 
013310 

01 332C 
013330 
013340 
01 3350 
013360 
013370 
013350 
013300 
0 1 3400 
013410 
013420 
013430 
0134UC 
013450 
013460 
01347C 
C134''C 
013490 
013500 
013510 
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** V ( J, K ) =V( J, K I 4U ( J,L »*VVf L,K ) 
RETCRN 
ENO 



SUB.’CUTIN'E CELVi(LUfSVI liVEtV) 

DIMENSION' UUI6,6)-,VV(6,6I,V(6,6).SVI (f.l I.SVEId.l) 
DO 9C M=),6 
1 F ( 'j V I ( N , II ) 1 , 90 , 1 
1 M=M+] 

DC 91 J= 1 ,6 

91 VV( J,M)=UU( J,M 
90 CONIINur 

L = 0 

DC 9? N=l,f> 

IFISVCIN, 1 ) )9r,r,92 
2 L=L4l 

or 93 K=1 ,3 
93 VIL,K)=VV(N,K ) 

92 CONI INUE 
RETURN 
END 

END 



0 13520 
013530 
0 1351*0 
013550 
013560 
013570 
0135 EC 
013590 
0 1 3600 
013610 
01362C 
013630 
013640 
013650 
013660 
013670 
013680 
013690 
013700 
013710 
013720 
013730 
013740 
013750 
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appendix e 



ACCURACY OF I'lSTHOD 



E.l General 

Appendix E contains the vibration analyses of the 
several typical piping configurations to which reference 
is made in Chapter 3* The analyses are grouped according 
to the purpose for which they were performed. The firsi 
group vjas performed to establish the accuracy of the method 
embodied in the program VIPIPE, the second and third to 
establish its integrity. 

E.2 Method Accuracy Analyses 

E.21 Systems Analyzed (straight sections) 



System 


Systems 


Parameters 




1 2. 


1 3 . 


end conditions 


1 t 

fixed-fixed 


^fixed-free 


fixed-cr op'p e'd 


length (in.) 


200 


200 


200 


diameter (in.) 


2.0 


2.0 


2.0 


wall thickness (in.) 


.125 


.125 


.125 


density (Ib./cu.ft.) 


1172 


1172 


1172 


shear modulus (psl) 


30 * 10^ 


30 * 10^ 


30 *10^ 


elastic modulus (psi) 


12* 106 


12- 10^ 


12. 10^ 



E-1 



^o 22 Goiiparlson of Natural F.reouencle! 






3 .C System 1 



Mode 

number 


frequency (rado/seco) 


Difference {fo) 


VIPIPE 


Comparison Values 


1 


48 o 5330086 ? 


48,533126 


0 00024 


2 


133 c ?8346278 


133.78345 


oOOOOO 


3 


262 o 26843097 


262 o 26886 


0OOOI6 


4 


^33.54478122 


433 054393 


0OOOI9 


5 


647062290622 


647.63991 


0OO26 



bo System 2 



Mode 

number 


frequency (rado/seco) 


Difference (^) 


VIPIPE 


Comparison Values 


1 


7.62356485 


7.6270709 


.04595 


2 


47.79657598 


47.797761 


0 00247 


3 


133.8363688 


133.83635 


oOOOOl 


4 


262c26528107 


262 026544 


.00006 


5 


433.54435149 


433.54415 


.00004 



Co System 3 



Mode 

number 


frequency (rad./sec.) 


Difference 


VIPIPE 


Comparison Values 


1 


33.44525745 


33.445600 


.00102 


2 


108.38602303 


108.38599 


.00003 


3 


226,1387715 


226.13868 


.00000 


4 


386,71069206 


386.71057 


,00003 


5 


590.10235871 


590.10164 


.00012 



^Sources of comparison frequencies for problems 1 
through 3 are given in Eo23o 



E,23 Source of Conparlson Frequencies for oyste.us 1 through 3 



The governing fourth order equation of a horaogeneous 
single cooiponent straight section system 



y,lv = UCd2 



vjas solved to give 

y = a*coso<rx + b*coshocx + c*sinpCx + d*slnhccx 
Upon applying boundary conditions for problem 1 
(y=y^ =0 at x=0, and x=L) and simplifying It was found 
that the eigenvalues of the system must satisfy 

cos oc CO shoe =1 (E.23.I) 

where 




Upon applying the boundary conditions of problem 2 
(y=y'=0 at x=0,y“=y”*=o at x=L) It was found that the 
eigenvalues of the system must satisfy 

cosoCcoshoc=-l (E»23"2) 

Upon applying the boundary conditions of problem 3 
(y=y'=0 at x=0 ,3E=y't=0 at x=L) It was found that the 
eigenvalues of the system must satisfy, 

sinoccoshoc -cosocslnhoc=0 (E»23-3) 

The roots of equations E,23-l, E,23-2 and E.23-3, and 
from them the first five natural frequencies of the systems, 
vjere obtained through the use of a digital computer. 



iLcZk- System Analyzed (section of constant curvature) 



System para meters % 



radius of curvature (ino) 


50 


included angle of arc 


270° 


diameter (in^) 


2 


wall thickness (in®) 


.125 


elastic modulus 


30 » 10^ psi 


shear modulus 


12 » 10^ psi 


density (lbs c/cu^-ft o ) 


556 




^ ,, end conditions? fixed-^fixed 

Results: 

2 

in-plane vibrations 



Mode 

number 


f*]?GQii6ncy (ro-d.o/s6Co ) 


Difference (^) 


VIPIPE 


Comparison Values 


1 


73o32973072 


73o6 


,23 


out-of-plane vibrations 


Mode 


frequency (rad ./sec.) 






VIPIPE 


Comparison Values-*- 


Difference {%) 


number 








1 


38.87735138 


39,8 


2,3 


Ec3 


Method Integrity Analyses 





The first group of analyses considered here constitutes 
a mirror Image check on system integrity o Four configura- 
tions were considered, and several mathematical models 
were usedo The second group consists of two analyses of 



This value was obtained using a curve in reference 2 
which 'was constructed through the use of a modified Rayleigh 
methodo In the words of the reference^ " o o =. laboratory tests 
have shown the curves to be substantially correct 

^ Comp arison values obtained by employing the equation 
co=F(o^-^EJ^ after obtaining F from curves in ref<, 2o 




the same single branch system v;hereln first one member then 
another is considered to be the branch. 

E. 3 I Mirror Imap:e Analyses 
a. System 1 

system 

mirror image 

System parameters: See system 2 E.21 

Mathematical model: Distributed mass with the 

effects of shear deflection and rotational Inertia considered. 




Results : 



Mode 

number 


frequency ( rad ./sec.) 


Difference (^) 


System 


Mirror Image 


1 


7.622895556 


7.62289556 


.00000 


2 


47.7668I87I 


47.7668I872 


,00000 


3 


133.63614681 


133.63614691 


.00000 


4 


261.54431793 


261.54431024 


.00000 


5 


431.64844506 


431.64917860 


.00016 



b. System 2 



System parameters: 



system 

77f7 

mirror image 
See system 3j E.21 
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Mathematical model: Distributed mass with the 



effects of shear deflection and rotational inertia neglected® 



Results : 



Mode 

number 


frequency (rad. /sec,) 


Difference {%) 


System 


Mirror Image 


1 


33 .W 25745 


33*445257^5 


.00000 


2 


108.38602303 


108,38602303 


,00000 


3 


226.13877715 


226,13877715 


,00000 


4 


386.71069206 


386.71059481 


.00002 


5 


590.10235871 


590.08973387 


,00211 






Component 


Parameter 


1 


2 


3 


4 


5 


6 


7 


diameter (in.) 








H 








2.0 


2,0 


2,0 


A 


2.0 


2.0 


2.0 




wall thickness (in.) 


.125 


,125 


,125 


N 


,125 


*125 


,125 


length (in.) 


130. 


1,40 


6,0 


G 


60,21 


1,50 


6,28 


radius of curv(in.) 


cc? 


4.0 


cxO 


E 


30. 


00 


4,0 


incl, angle of arc 


0 


0 

0 


0 


R 


115*° 


0 


90° 



E -6 





ntfck 



Intensive properties of all co:nponsnts except 

hangers : 

density: 556 lbs,/cu.-ft. 

shear modulus: 12 • 10^ psi 

elastic modulus: 30* 10^ PSi 

Hanger spring rates: 

CLX = 1000 lbs. /in. 

CLI = CLZ = 1000 lbs, /in. 

CTX = CTI = CTZ = 1000 in. -lb. /rad. 

End conditions: 

Both ends of system fixed. 

Mathematical model: Distributed mass with the 

effects of shear deflection and rotational inertia considered. 



Results: 

in-plane vibrations 



Mode 


frequency (rad. /sec.) 




number 


System 


Mirror Image a. 


Difference (^) 


1 


116,91660080 


116.91660084 


.00000 


2 


305.21721131 


305.21720649 


.00000 


3 


567.56373058 


567.56397808 


.00005 


4 


845.99432421 


846.17624950 


.02150 


5 


— 


1135.09717500 
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out-of-plane vibrations 



node 

number 


frequency ( rad » /sec,) 


Difference (^) 


System 


rlirror Image a. 


1 


114,21818694 


114 c 21818694 


,00000 


2 


243,92831220 


243.92831317 


.00000 


3 


321,40891123 


321,40891770 


.00000 


4 


591.40446121 


591.40449774 


,00000 


5 


872,23102583 


872,24253187 


,00131 


6 


I1380OI776773 


II380I9615015 


,01567 


7 


1571«78173050 







in-plane vibrations 



Mode 


frequency (rad, /sec,) 




number 


System 


Mirror Image b. 


Difference (/o) 


1 


116,91660080 


116,91660080 


.00000 


2 


305.21721131 


305.21721131 


oOOOOO 


3 


567.56373058 


567.56373058 


,00000 


4 


845.99432421 


845.9943421 


,00000 




out-of- 


plane vibrations 




Mode 


frequency (rad, /sec,) 




number 


System 


Mirror image b. 


Difference (/a) 


1 


114,21818694 


114.21818694 


,00000 


2 


243.92831220 


243.92831220 


,00000 


3 


321 ,40891123 


321.4089II23 


,00000 


4 


591.40446121 


591.40446121 


,00000 


5 


872.23102583 


872,23102583 


,00000 


6 


1138.01776773 


1138,01776773 


,00000 


7 


157I.78I73050 


1571.78173050 


oOOOOO 





system (component 
numbers shown) 



^jgr^indlcates pipe 
hanger 



mirror image 





Component 


Parameter 


1 


2 


3 


4 


5 


6 


7 


diameter (in.) 


2.0 


2.0 


2.0 


2.0 


H 


2.0 


2.0 


wall thickness (in.) 


.125 


.125 


.125 


.125 


A 


.125 


.125 


length (in.) 


400. 


39.2? 


100. 


39.27 


N 


100. 


100. 


radius of curv.(in.) 


oo 


50. 


csO 


50 . 


G 


exO 


CxO 


incl. angle of arc 


0° 


^50 


0<^ 


450 




0° 


0° 





Component 


Parameter 


8 


9 


10 


11 


12 


13 


diameter (in.) 


2.0 


H 


2.0 


2.0 


2.0 


2.0 






A 










xvall thickness (in.) 


.125 


N 


.125 


.125 


,125 


.125 


length (in.) 


39.27 


G 


39.27 


100. 


39.27 


100. 


radius of curv.(in.) 


50 . 


E 


50 . 




50 . 




incl. angle of arc 


450 


R 


45O 


0° 


450 


QO 
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Incensive orooertles of all co'nponents except 

hangers : 

densitys 460 IbSo/cuo-fto 
shear modulus: 12 » 10^ psi 

elastic modulus 30 - 10^ psi 

Hanger spring rates: 

CLX = CLI = CLZ = 1000, IbSo/lno 
CTX = CTY = CTZ = 1000, in. --lb , /rad. 

End conditions: all main member and branch ends 

fixed. 

Mathematical model: Distributed mass with the 

effects of shear deflection and rotational inertia considered. 



Results : 

in-plane vibrations 



Mode 

number 


frequency (rad/sec,) 


Difference {%) 


System 


Mirror Image 


1 


18,55192182 


18,55189208 


0OOOI6 


2 


5I02521I856 


51,25212736 


,00001 


'i 








j 








4 




107,19^59215 
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out-of-plane vibrations 



node 

number 


frequency (rad. /sec.) 


Difference(/a) 


System 


Mirror image 


1 


14,31754957 


14,31754879 


,00000 


2 


34.45495449 


34.45494654 


,00002 


3 


41.07914530 


41.07913499 


.00002 


4 


73.06452245 


73.06429798 


.00030 


5 


115.86572509 


115,86505834 


,00057 






Component 


Parameter 


1 


2 


3 


4 


diameter (in.) 


2.0 


2,0 


2,0 


2,0 


wall thiclmess (in,) 


,125 


,125 


.125 


.125 


length (in.) 


100 . 


100 . 


,01 


100 . 


radius of curv,(in.) 


oO 


c =-0 


.01 


oO 


Incl. angle of arc 


QO 


0 ° 


45° 


0 ° 
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Intensive properties of components: 



density: 55^ lbs e/cuo-ft « 

shear modulus: 12<= 10^ psl 

elastic modulus: 30“ 10^ psi 

End conditions: Both ends of main member, and 

end of branch, fixed. 

Mathematical model: Lumped mass with the effects 

of shear deflection and rotational inertia neglected. 



Results : 

out-of“plane vibrations 



Mode 

number 


frequency ( rad . /s ec . ) 


Difference {%) 


System 1 


System 2 


1 


49,64604938 


49.64604961 


,00000 


2 


135.52684605 


135.52685066 


.00000 


3 


269,83567590 


269.83569992 


.00000 


4 


530.72012316 


530.75718901 


.00698 


5 


675.51662199 


675.57044008 


.00796 


6 


901.52786818 


901,76873827 


.02671 
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In-plane vibrations: 



Mode 

number 


frequency (rado/sec.) 


Difference (^) 


System 1 


System 2 


1 


132.89829688 


132.89836401 


.00005 


2 


190.61015386 


190.61009901 


.00002 


3 


430.03806330 


430.03948554 


.00033 


4 


512.55370411 


512,54463160 


.00177 


5 


893.66594383 


893.10993095 


.06221 


6 


962.07559083 


963.00709091 


.09682 
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iiPPi:;NDIX P 



TRAxspi:,R Matrix source RajParexce 



F.l General 

Fourteen of the transfer matrices used in VIPIPS 
were drawn from the catalog of transfer matrices, pages 
3?6-410, in Pestel's "Matrix Methods in Elastomechanlcs" 

Hill . Some of the transfer matrices constructed in VIPIPE 
are a composite of two separate transfer matrices from[l]. 
These composite matrices are constructed as shown in 
Figs. F.1-1, and F.1-2, The former shov/s the format of 
case 1 composite transfer matrices, the latter those of Case 2, 
In both 



^1,1 I 



^ 2,1 : 



° 1,2 



I 0 



2,2 



B 

null 

submatric^ 



Fig. F.1-1 In-plane composite 
transfer matrix format 



Fig. F.1-2 Out-of-plane 
composite transfer n.atrix 
format 



cases the a matrix is of fourth order, the B .natrix of second 
order. In Case 1 composite transfer matrices the elements 
of the B matrix occupy the positions shown (Fig. F.1-1). 



F-1 




I 



In the follovjing source listing, v^hen tvjo matrices are 
indicated for one subroutine, the first is the A matrix, 
(All component transfer matrices are constructed in sub- 
routines 0 ) 

P,2 In-Plane Vibration Transfer Matrices 



Transfer matrix Subroutine Source 



straight section with distri- 
buted mass, shear deflection 
and rotational Inertia 


DISTM 


C-3a 

C-la 


straight section, elastic 
massless field with shear 
deflection 


SFIELD 


C-3c 

C-lc 


curved section, elastic mass- 
less field without shear 
deflection 


C FI ELD 


C-7c 


point mass vjith rotational 
inertia 


POINT 


c-3f 

C-lb 


rigid straight section with 
rotational inertia, distributed 
mass 


RIGID 


C-3d 

C-lb 


stiff corner (elbow), massless 
field 


STIFCO 


C-7f 


Out-of-Plane Vibration Transfer Matrices 




Transfer matrix Subroutine 


Source 


straight section with distri- 
buted mass, shear deflection 
and rotational inertia 


DIDTMO 


C-3a 

C-2a 


straight section, elastic mass- 
less field with shear deflection 


3FIEI0 


C-3c 

C-2c 


curved section, elastic mass- 
less field with shear deflection 


CFIELO 


C-8c 


point mass with rotational 
inertia 


POINO 


C-3f 

C-2b 


rigid straight section with 
rotational inertia, distributed 
mass 


RIGIO 


C-3d 

C-2b 


stiff corner (elbow), massless 


STIFOO 


C-8f 



field F‘2 



AFPEIMDIX G 



oAi'lP LE F ROBLiiiM 
G,1 General Remarks 

Appendix G. is Intended to supplement apoendix B, 

Included in appendix G. is an illustration of the data 
deck for a typical multiple branch system with hangers, 
and an example of the three output formats available in 
VIFIPE. 

G.2 The Data Deck 

The data deck as shown calls for determination of both 
Case 1, and Case 2 natural frequencies (IoP=0). If only 
Case 1 frequencies were sought, card 4 would be removed 
from the deck, and lOP^ (circled field on the first data 
card) would be set equal to 1. Similarly, if only Case 2 
frequencies were sought, cards 2 and 3 would be removed, 
and lOP set equal to 2. 

In the sample problem configuration all boundary con- 
ditions may be specified by code (Bo6l), but in order to 
Illustrate the method of programming boimdary conditions 
for which no code is provided, the Case 1 boundary condition 
state vector for branch 1 is read-ln directly o Note that 
the second field of card 2 is blank (blanks are read as 
zeros by the computer) and that card 3 contains the particular 

^See appen, B<,, section B,71« 



B t'dnch 




, . 1 . . ; . 



Fig. Gr.2-1 Multiple Bx-anth S^s 



44.0. 30.00E'06nr2T00E 04 

“27 -.125 - loor 

“27 7125 -50. ■ 45. 

a. .125 100. 






1 



* 0 !► 



boundary condition state vector o 

The data deck used in Fle» G«2-2 pertains to system 4 
of section Eo31> apnendix E<, The system is illustrated in 
r’i,^. G. 2 -I 0 
G .3 Output Options 

Examples of the three output formats available in 
VTPIPE, through specification of the value of the parameter 
GJi:iBC (B. 7 I), are included on pages G-4, and G-5«. The 
system analyzed was system 3 of section E. 31 , appendix E. 
and is illustrated there* The problem numbers refer to 
the order in which the data decks are "stacked" (B 08 ). 

The heading ... "PROGBiiM VIPIPE" ... etc. is printed once 
per run, e.g., if 4 data decks were "stacked" behind the 
program deck, and each asked for Case 1 and Case 2 analyses 
the heading would precede the print out of the results of 
the Case 1 analysis of system 1 only. 
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M 0 G 14 K vlf>IPE C. E. FINK NFA-I ' FEEHUARV 1 «N 

PETMOO of transfer MATR^CE^ FRESLENCIES Df A PLANAR PIPING SrSTtH ARE OITERPINEO BY AN irERAtlVE PROCECURE USING THE 



IN PLANE PCOE FREOLENCIES IRAOIANS PER SECCNCI 
MODE FREGLENCY 

1 11 «. 91 « 60 a 8 k 

2 }os.2iTaca«2 

3 '• SET.SEREsno 

4 a 4 «. 1404 S 6 S 1 

PIOaLEN TERMINATED CUE TO SIGNIFICANT FIGURE LIMITATION OF CCMPUTER. 

efle^iion'^ano roIar*’’"*"^- 



E effects op shear OEFLE 



i»?A“ilR|’'?feSlf8E«8. 



Fig. G.3-1 GDABC=Q 



aWMAM IFfPIPC C. E. 

„ . MODAL FPEDUENCIES 3 P a planar pipin 

NETNOO pr transper matrices. 



PROBLEM t 

IN plane mode FRFBUENCIES IRAOIANS PER lECONO) 
NODE PREOUENCV 

1 nt.aiMaoso 

S . SC5.21T2I0T3 

a SAT.SASaSITA 

4 eua.osioTiaa 

a 1l28.3T4S4a7» 

PROBLEM TERMINATED DUE TO SIGNIFICANT FIGURE LINITATICN OF CCMPUTER. 

THE EPFECTS 6 f SHEAR DEFlIIt liS^ISo ROTARY^I !!iR?llA‘‘C iRl'^ciSNa ISeReS. 

fNP 0 |* 0 ATA 

COMPONENT DIAMETER 



rUCRNESS tENGTH 



fS8»lu8^ 






mil 

9 o: 8 s 



MgfeSfJi' 

•yiii 



MOOUL&S I 



|O^NOARY CONOljtONS ^ 
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